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ABSTRACT 


The construction of an apparatus for the condensation 
performance testing of a horizontal bundle of four tubes with 
various refrigerants was completed. The apparatus was 
instrumented, and data reduction software was developed to 
provide bundle and single tube condensation data. 

Two tube bundles were tested, smooth copper tubes and low 
integral-fin copper-nickel tubes, with two refrigerants, R-114 
and R-113. An enhancement ratio of about 2.0 for the overall 
heat transfer coefficient was demonstrated for the finned 
tubes over the smooth tubes. 

Internal contamination, possibly due to a breakdown of the 
refrigerant molecules when subjected to high temperatures in 
the boiling chamber, inhibited further meaningful data 
Gorrection. Recommendations for improvement of the test 


apparatus are made. 
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A. BACKGROUND 

With the advent of more complex shipboard weapons/combat 
systems and increasing fuel costs, the United States Navy has 
increasingly recognized a need for an energy-efficient, light- 
weight, and high-capacity air-conditioning system. Such an 
advanced system has been proposed and is currently in its test 
and evaluation phase for the DDG-51 program [{Ref. 1]. 

Among the substantial differences in the advanced air- 
conditioning plant over extant systems is the use of titanium 
finned tubes in the condenser, in place of the smooth copper- 
nickel alloyed tubes previously used. The use of titanium in 
the air-conditioner's condenser presents two significant 
advantages. First, since the United States Navy uses ambient 
seawater as a heat sink, the corrosion resistance of titanium 
should measurably improve system integrity and reliability. 
Secondly, the advanced air-conditioning system expects to 
realize a weight savings per unit of 2286 kg (5040 lbs) and 
the major portion of this weight savings will come from the 
use of the lighter-weight titanium condenser tubes [Ref. 1]. 

Since the early 1970's the United States Navy has used R- 
114 as its primary refrigerant. The advantages of R-114 over 
the more-widely-used commercial refrigerants in a naval 


application derive from its low toxicity, temperature 


stability, stability when in contact with moisture, and its 
applicability to lower-pressure systems [Ref. 1]. R-114’* will 
be the refrigerant utilized in the advanced air-conditioning 
milanies 

As no substantial data base on the heat-transfer capacity 
exists for a system utilizing titanium condenser tube bundles 
configured as proposed in the advanced air conditioning plant 
and consequently, no way exists to measure this design's 
performance against competing designs, the need arose to 
develop a test apparatus to provide this data base in order to 
accurately predict future system performance. The design and 
subsequent construction of this test apparatus at the Naval 
Postgraduate School began in 1987 and is documented in the 
Theses of LCDR David S. Zebrowski [Ref. 2] and of LCDR Thomas 
J. Murphy [Ref. 3]. The test apparatus was also intended to 
serve as a test platform for advanced boiling surfaces 
proposed for use in the advanced air-conditioning plant. The 
incorporation of the advanced boiling surface test platform 
and a condenser tube test platform into a single test 
apparatus allows the widest latitude in examining the effects 
of various heat inputs, coolant flow rates, and various levels 
of refrigerant contaminants on individual components’) and 


overall system performance. Zebrowski and Murphy constructed 


‘New refrigerants may have to be used in the future due to 
the "ozone problem." 


the basic apparatus, and Murphy operated the system for 


preliminary evaporator measurements. 


Bs OBJECTIVES 
The major objectives of this thesis were: 


1. Refine the design and complete the fabrication of the 
apparatus for the testing of condenser tube bundles 
utilizing R-114 as the working fluid. 


2. Develop and instrument a data-acquisition system with 
the associated software to provide a data base on the 
heat-transfer capacity of the condenser section of the 
test apparatus. 


3. Validate the test apparatus by comparing condenser 
performance against extant data bases derived from 
conventional condenser tube bundles. 


II. LITERATURE SURVEY 


A. GENERAL OBSERVATIONS 

Condensation is the phase transformation process in which 
a vapor is transformed to liquid by removal of latent heat. 
The promotion of condensation in heat exchangers is used 
extensively in applications for propulsion engineering and 
air-conditioning/refrigeration cycles. Due to its importance 
in the aforementioned fields, considerable research has been 
directed at the factors influencing the process and ways in 
which the process can be enhanced. Such factors as various 
modes of condensation, surface orientation to vapor flow, the 
shear forces exerted on condensate film, various external and 
internal surface enhancements, turbulent effects due to 
cascading condensate flow from another surface, and the 
effects of fluid properties on the process are detailed 
extensively in various reviews. This thesis deals exclusively 
with film condensation of refrigerants on small horizontal 
tube bundles in a quiescent vapor where tubes with external 


fins are compared to smooth tubes. 


B. SINGLE TUBE INVESTIGATIONS 
1. Smooth Horizontal Tube Studies 
The first comprehensive condensation model was 
developed by Nusselt in 1916 [Ref. 4] based on the assumption 


that a quiescent vapor at saturation temperature coming into 


contact with a wall surface below saturation temperature would 
condense and form a continuous film of condensate growing in 
thickness as the film flowed off the surface under the 
influence of gravity. The condensate film at the vapor-liquid 
interface would be at the vapor saturation temperature with a 
temperature gradient in the film down to the wall surface 
temperature. No radiation or convection would take place at 
the liquid-vapor interface as both liquid and vapor are at the 
Same temperature, however the amount of vapor condensing 
corresponds to the quantity of heat flowing through the film 
by thermal conduction. In the case of a horizontal tube, he 
assumed that condensate flows in a laminar manner around the 
Sides of the tube and off the tube in a continuous sheet. The 
heat transfer coefficient would be maximum at the top center 
of the tube decreasing gradually around the surface of the 
tube, as the inherent thermal resistance in the film grows 
with the thickness of the film, and eventually goes to zero at 
tube bottom. Nusselt's model is limited by disregarding 
surface tension forces that tend to hold up condensate at the 
tube bottom until overcome by gravity, resulting in the 
production of condensate droplets rather than a continuous 
sheet. Nusselt's expression for the average heat transfer 
coefficient for a single tube subjected to a constant heat 


flux is given by: 


_ 3. 2 ae ye (2a 
a, = -6551\--p_e-g-Ah /n_-D, "0" 


Although constrained by the aforementioned limitations, 
Nusselt's model remains the conservative bench-mark against 
which all other models are compared. 

2. Exterior Surface Enhanced Tubes 

Research into ways to improve condensation performance 
in condensers is motivated by the idea that any force that 
acts to thin the condensate film promotes an enhancement of 
heat transfer by minimizing the resistance to heat flow. 
General enhancement techniques include the fabrication of 
surfaces that promote dropwise condensation, wrapping exterior 
surfaces with wire, installation of porous drainage strips, 
and the fabrication of finned tubes with fins of various 
geometries and various spacings along the tube. This thesis 
deals with a low integral-fin tube and its heat transfer 
enhancement over a smooth tube. 

In considering a low integral-fin tube during 
condensation, there exists two distinguishable regions on the 
circumference of the tube; a flooded region and an un-flooded 
region. The flooded portion of the tube defines’ the 
condensate retention angle of the tube with respect to the 
tube circumference. The smaller the condensate retention 
angle of the tube, the larger the heat transfer capacity of 
the tube. The fin can be divided into three regions; the fin 
tips, the sides of the fins, and the fin root area. The 
majority of heat transfer takes place at the fin tips. 


Surface tension forces pull the condensate from the tips down 


the fin sides into the flooded root area, where gravity drains 
the condensate to the bottom of the tube. The amount of 
condensate retained along the circumference of the tube is 
dependent upon the ratio of surface tension forces to gravity 
forces. 

Beatty and Katz [Ref. 5], in 1948, developed a 
comprehensive model for the prediction of the heat transfer 
coefficient of finned tubes based on their experimental 
results for various test fluids (including R-22) and finned 
tubes of various metallic compositions, and various fin 
geometries. Their model assumes gravity-dominated flow and 
neglects surface tension effects completely. Subsequent 
experimental results from other sources indicate that the 
Beatty and Katz model overpredicts the heat transfer 
coefficient as surface tension increases or as fin density 
increases. Nevertheless, the Beatty and Katz model conformed 
to their reported overall heat transfer enhancement of up to 
2.3 for R-22 on finned tubes compared to smooth tube data. The 


Beatty and Katz model is given by: 


~ Saee2 ; 1/4.,1,1/4 72 OS 
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A three-region theoretical model, based on their experimental 
results with three rectangular finned tubes, was presented by 
Karkhu and Borovkov [Ref. 6], in 1971, including surface 
tension forces. Measured heat-transfer coefficients 
demonstrated a 50 to 100 percent increase in vapor-side heat 
transfer coefficients for steam and R-113 condensing on finned 
tubes compared to smooth tubes. Unfortunately, they did not 
report enhancements separately for the two fluids. In 1980, 
results reported by Carnavos [{Ref. 7] condensing refrigerants 
on various finned tubes demonstrated an enhancement of up to 
400 percent in the heat transfer coefficient over smooth tube 
results. Work done by Sauer and Williams [{Ref. 8], in 1982, 
on the condensing performance of finned tubes with oil- 
contaminated R-113 demonstrated a serious degradation of 
performance when the surface tension to density ratio was 
large. The conclusion formed was that the higher-surface- 
tension oil remained in the fin gaps rendering the finned 
surface ineffective. Results reported by Honda et al. [Ref. 
9], in 1983, condensing R-113 on finned tubes of various 
geometries showed improvements of 900 percent for the vapor- 
side heat transfer coefficients at a constant vapor-to-tube- 


wall temperature difference. The results reported by Kabov 


[Ref. 10] in 1984 with refrigerants R-12 and R-21, indicated 
that the bulk of the latent heat was in fact removed on the 
lateral fin surfaces and that an optimum fin height and 
Spacing was dependent upon the ratio of surface tension to 
gravity forces. The results published by Masuda and Rose 
(Ref. 11] in 1985 in experiments with R-113 condensing on low 
integral-fin tubes, confirmed that the overall heat transfer 
coefficient increases, in general, with decreasing fin 
spacing. Their results showed a 600 percent increase in 
enhancement over smooth tube performance. More recent results 
by Marto et al. [Ref. 12] with R-113, demonstrated a 700 
percent enhancement in performance over a smooth tube and gave 
an optimum fin spacing of 0.5 mm for that fluid. Work 
performed and published in the same time frame by Sukhatme et 
al. [Ref. 13] with R-11 on conventional integral-fin tubes and 
special pyramid-shaped fin tubes, reported enhancement ratios 
of 5 to 7 for the low integral fin tubes and 10.3 to 12.3 for 


the pyramid-shaped fin tubes. 


C. TUBE BUNDLE INVESTIGATIONS 
1. Smooth Tube Bundles 
In smooth tube bundles, two conflicting factors play 
a role in determining bundle performance. First, the 
condensate flowing from the tubes above a given tube in a 
bundle tends to thicken the condensate film on that tube, 
hence increasing the resistance to heat transfer. This effect 


is known as the condensate inundation effect. Secondly, 


droplets from other tubes striking the film surface on a tube 
with a velocity provided by gravity or vapor flow, can create 
ripples or waves in the condensate film imparting a turbulence 
within the condensate film that produces an enhancement of 
heat transfer performance. 

In 1949, Jakob [Ref. 14] elaborating on Nusselt's 
model, predicted that the vapor-side heat transfer coefficient 
for a tube in a bundle, compared to the first tube in the 
bundle, was a function of that tube's relative position in the 
bundle. This model was based upon the assumption that a 
continuous laminar sheet of condensate flowing off the tube 
directly above the tube considered and striking the top of 
this tube further thickened the condensate film on this tube. 


Jakob's model is given by: 


all =wil/4 (2.4) 
1 

In 1958, Kern [Ref. 15] proposed a model based on the 
assumption that discrete droplets or jets of fluid from other 
tubes caused ripples in the condensate film diminishing the 
inundation effect. Kern's model is less conservative than 
Nusselt's and, in many cases, remains the applied industrial 
design standard. Kern's model is given by: 


a 
N = y 76 (225) 
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Work by Chen [Ref. 16], published in 1961, proposed a model 
that considered the momentum gain of falling condensate as 
well as the condensation of vapor on sub-cooled condensate 
droplets or sheets. Chen's model which is essentially 
Nusselt's model times a factor that incorporates the phase 


change number is given by: 


“NL yr l/4 


N [l + 2-Ph- (N-1) ] (26) 
a 


where the Phase Change number is given by: 


C AT 
Ph = — 


Ah. 





Equation (2.6) is valid provided, 


Pe aN ele 210) 


A model proposed by Eissenberg [Ref. 17] in 1972, assumes that 
condensate does not always drain in a vertical direction but 
may be diverted sideways due to local vapor flow conditions. 
In this case, the condensate strikes subsequent tubes on their 
Sides rather than on their tops, minimizing inundation effects 
on the condensate film in the top portion of the tube. 


Eissenberg's model is given by: 


| Q | 
Zi 


om 


= 60 + .42-(n7 174) (a7) 
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2. Enhanced Tube Bundles 

Katz and Geist [Ref. 18] in 1948, studying six fin 
tubes in a vertical row, using R-12 among other working 
fluids, found that the effect of condensate inundation was 
over-predicted by Jakob's model and proposed that the exponent 
in Equation (2.4) be changed to 0.06 for finned tubes. 

A theoretical model was proposed by Honda et al. [Ref. 
19} in work published in 1987, for finned tube bundles that 
showed good agreement, within 5 to 7%, for data taken 
condensing acetone and R-12 on finned tubes. Their model 
relies upon solving a set of algebraic equations describing 
the vapor-to-coolant conjugate heat transfer problem. A 
compendium of the aforementioned models, both single tube and 
bundles, is provided by Marto [Ref. 20]. 

As of the date of this writing, the author is unaware 
of any comprehensive model, substantiated by experimental 
results, that accurately predicts heat transfer performance 
for enhanced tube bundles with varying fin geometries and 


varying pitch. 
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III. EXPERIMENTAL APPARATUS 


A. CONDENSER/BOILER TUBE BANK TEST PLATFORM 
1. Overview 

The composite test platform, with associated support 
systems, is depicted schematically in Figure 3.1. The boiler/ 
condenser unit was fabricated from 6.35 mm thick, rolled 
stainless steel plates designed to withstand an absolute 
pressure of 308 kPa. The top cylinder which serves as the 
condenser for the apparatus is 1.30 m in length with an 
external diameter of 0.61 m. The effective condensing length 
is 1.22 m. The condensing chamber is capped on either end 
with circular stainless steel plates of 0.71 m diameter and 
connected to the chamber with a flange and stud assembly. 
System integrity is maintained with a 3.20 mm thick rubber 
gasket. The end caps support the nylon block condenser tube 
mounts, stainless steel backing plates, auxiliary condenser 
coils (coolant entrance side), and mixing chambers (coolant 
exit side). The condenser chamber has five 12.7 mm thick 
Pyrex glass view ports backed with 12.7 mm thick Plexiglas. 
The view ports are 127 mm in diameter and are located axially 
along the main chamber, angled so as to provide top to bottom 
views of the test tubes at various locations along the 


effective condensing length. The condenser’ chamber and 
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ancillary equipment are depicted graphically and in 
photographs in Figures 3.2 through 3.5. 

Internal to the condenser chamber, a stainless steel 
shroud is fitted over the effective condensing length of the 
instrumented test tubes (see Figure 3.6). The purpose of the 
shroud is to channel refrigerant vapor along the inside 
circumference of the shell, collecting vapor at the top, and 
then forcing vapor through the vertical axis of the 
instrumented tube bank into the shroud's well, where 
condensate is collected and the remaining vapor is condensed 
by an auxiliary condenser. The auxiliary condenser is 
composed of five helically wound copper tubes of 9.53 mm 
diameter suspended inside the shroud well by cantilevered rods 
welded to the entrance end cap. The stainless steel shroud 
was fabricated with a glass panel serving as one side of its 
stem to permit viewing of the instrumented tubes through the 
view ports. 

The condenser chamber 1s attached to the boiler 
chamber by a rolled stainless-steel cylinder nominally 280 mm 
in diameter and 203 mm in length, located mid-way along the 
condenser chamber's length allowing condensate to drain by 
gravity. A collar dam was fitted at the connecting flange, 
for the initial purpose of providing a condensate drain point 
to send condensate to the auxiliary storage unit. 
Subsequently, a faster, more-efficient method of transferring 


refrigerant was devised that allowed refrigerant vapor to be 
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Figure 3.6 Cross-section of Condensing Chamber Schematic 
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sent and condensed in the auxiliary storage unit. A 9.53 mm 
diameter copper tubing was provided for the vapor flow from 
the top of the condenser chamber to the storage tank. A 
similar tubing was provided between the storage tank and the 
boiling chamber to allow liquid flow via gravity. 

The boiling chamber is connected to the condenser unit 
via the interconnecting riser as previously discussed. The 
boiling chamber was also fabricated from rolled stainless- 
steel and is cylindrical with an outside diameter of 0.61 m 
and a length of 0.279 m. The boiling chamber is fitted with 
two Pyrex view ports, strengthened by Plexiglas plates, for 
observation during operation. The boiling unit is comprised 
of three groups of tubes (see Figure 3.3). In the simulation 
tube bundle, there are five boiling tubes each nominally rated 
at 1.5 kW and located at the bottom of the boiler unit. In 
the auxiliary tube bundle, there are four boiling tubes each 
nominally rated at 4 kW and located two on either side of the 
instrumented tube bundle. The instrumented tube bundle is 
located in the lower center of the boiling unit and is 
comprised of 35 tubes, with ten active tubes rated at one kW 
each, five instrumented tubes, and 20 dummy tubes. The power 
provided to these tubes is controlled by three variac 
controllers located on a console adjacent to the test platform 
and fed by a 208 volt power supply. The variac controllers 
are graduated in one percent increments of maximum power for 


each tube bundle. The exact location of each tube bundle is 


20 


apparent in Figures 3.2 through 3.5. It is a unique feature 
of this test platform that both boiling and condensation 
phenomena in a closed loop system can be evaluated 
simultaneously. Details on the design and construction of the 
basic test platform are available from Zebrowski [Ref. 2] and 
Ma@rphy (Ref. 3]. Lightoff and securing procedures for the 
apparatus are listed in Appendix B. 
2. Ancillary Systems 

Component equipment that support the R-114 tube bundle 
test platform and are located external to the apparatus 
include: (1) an R-114 storage and transfer system, (2) a 
condenser cooling and flow control system, (3) a coolant 
(62.4% by weight mixture of ethylene glycol and water) sump, 
and (4) an eight-ton refrigeration unit. 

The R-114 storage and transfer system, as previously 
described, consists of a stainless-steel cylindrical tank 
0.350 m in diameter and 0.91 m in length located on a rack 
above the coolant sump. Transfer of the refrigerant is 
accomplished by boiling in the main boiling chamber with vapor 
being sent to the storage tank via 9.53 mm diameter copper 
tubing located in the top center of the condenser chamber. 
The vapor is condensed in the storage cylinder by means of a 
helical copper coil, suspended the length of the cylinder on 
a cantilevered bar, that is kept cooled by the water-ethylene 
glycol mixture. Liquid refrigerant can be returned by gravity 


to the test platform from the bottom of the storage tank 
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through a 9.53 mm diameter copper tubing to the boiling 
chamber. Faster transfer of R-114 liquid was also possible if 
the system pressure was at or below the atmospheric pressure. 
Notice that the storage tank experiences an absolute pressure 
Of about 2107 KPa: 

Coolant and flow control to the test apparatus and 
ancillary equipment is accomplished by two different flow path 
systems. Both flow systems are driven by two 0.5 HP constant- 
speed pumps that take suction on the main coolant sump. The 
coolant for the test condenser tube bank passes from the pump 
Gischarge through 76 mm diameter PVC piping to a Plexiglas 
header. At the header, flow is split and proceeds through 
15.9 mm diameter Tygon flexible tubing to a bank of 
rotameters. Flow can be controlled by throttling a gate valve 
located at the entrance to each flow meter. Coolant flow 
leaves the exit of each flow meter and proceeds to the 
instrumented condenser tube bank through flexible Tygon 
tubing. At the exit of the main condenser chamber, coolant 
leaves each tube and flows through flexible Tygon tubing to 
individual mixing chambers (Figure 3.7). From the mixing 
chambers, coolant passes through flexible Tygon tubing to a 
central Plexiglas header suspended above the test platforn, 
and the collective coolant is piped back to the main sump 
through 76 mm diameter PVC pipe. Auxiliary coolant flow 
leaves the pump discharge through PVC piping and proceeds to 


a two-way ball valve, where the flow is either sent to the 
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Figure 3.7 Mixing Chamber Schematic 
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R-114 storage tank and back to the main sump or is sent 
through PVC pipe to a large rotameter. From the auxiliary 
flowmeter, coolant enters a central header and is then split 
to flow to the five auxiliary condenser coils and back to the 
coolant sump. Flow through each auxiliary condenser coil is 
controlled by a valve located at the coil inlet. The coolant 
sump has a 1.81 cubic meter capacity and is constructed of 
12.7 mm thick sheets of Plexiglas. The coolant is 
approximately 62.4% ethylene glycol and 37.6% water (by 
weight). 

The coolant is chilled by an externally-located, 8 ton 
refrigeration system that continually re-circulates sump 
coolant with a 0.75 HP pump. The refrigeration system is 
capable of maintaining a sump temperature between -21 C and 
ambient temperature. 

3. Instrumentation 

The coolant temperature rise is measured by series- 
connected thermopiles having ten junctions on either end. 
These thermopiles were fabricated using type-T thermocouple 
wire. The ends of the thermopiles are inserted in stainless 
steel wells with copper-plugged tips located at the entrance 
to each individual condenser tube and at the exit of this 
tube's mixing chamber. Inlet coolant temperature is measured 
with type-T thermocouples also located in the entrance wells. 


Refrigerant liquid and vapor temperatures are measured with 
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type-T thermocouples inserted in stainless-steel wells at the 
jocations indicated in Figures 3.2 through 3.5. 

System pressure iS monitored through a calibrated 
pressure-vacuum gage valid over a range of 30 inches mercury 
(vacuum) to 30 pounds per square inch (gage pressure). 

4. Tube Bundle Data Acguisition and Reduction 

A Hewlett-Packard 9816A computer was used to control 
a Hewlett-Packard 3497A Automatic Data Acquisition Systen, 
which read the output of the thermopiles and thermocouples. 
Readings were made in millivolts and were converted to 
temperature readings in the data reduction program. The 
channels read by the data acquisition system are listed in 
Table 3.1. 

5. Tubes Tested 

Two sets of tubes were tested. The first was a smooth 
copper tube (inside diameter 13.26 mm, outside diameter 15.88 
mm) and the second was a low integral-fin copper-nickel tube 
(inside diameter 10.16 mm, root diameter 14.00 mm, outside 
diameter 15.88). These two sets of tubes were tested, in 


bundles and individually, with R-114 and R-113. 
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TABLES3s 34 


CHANNEL ASSIGNMENTS ON DATA ACQUISITION SYSTEM 


Channel Numbers Measurement 

O=Z Vapor Saturation Temperature 

bo Liquid Temperature 

5 Inlet Temperature coolant list Tube 
6 Inlet Temperature coolant 2nd Tube 
7 Inlet Temperature coolant 3rd Tube 
8 . Inlet Temperature coolant 4th Tube 
20 Thermopile ist Tube 

on Thermopile 2nd Tube 

22 Thermopile 3rd Tube 

Ze Thermopile 4th Tube 
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IV. SYSTEM OPERATION AND DATA REDUCTION 


A. TEST PLAN AND MODIFICATION OF DESIGN FOR CONDENSER 
TEST APPARATUS 


1. Ability to Change Working Fluids 
Although envisioned in the original apparatus design 
by Zebrowski [Ref. 2] and Murphy [Ref. 3], no physical system 
existed for the change-out and storage of working fluids. The 
system was designed as a general test platform to evaluate the 
performance of various refrigerants and steam during 
condensation on a variety of test tubes. 
The R-114 storage and transfer system as described in 
Chapter III, ancillary systems, was designed and built to 
allow for the storage of refrigerants that evaporate at 
atmospheric temperature and pressure. This additional 
capability not only conserves expensive refrigerants but 
facilitates speedy change of these working fluids during the 
evaluation of a particular tube with a variety of 
refrigerants. 
2. Tube Alignment and Reduction of Tubes in Test Bundle 
Visual inspection, during preliminary experimental 
runs, revealed that bowing of the tube bundle _ was 
detrimentally affecting condensation patterns. As a result of 
misalignment of the endcaps during fabrication, bending 


moments of varying magnitudes and directions resulted in 
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condensate flow striking subsequent tubes in the vertically- 
oriented bundle at different positions around the 
circumference of the tubes. Upon disassembly of the end 
plates of the apparatus and comparison with the proposed 
Grawings for the DDG-51 refrigeration condenser provided by 
the David Taylor Research Center, miscalculations in the tube 
pitch were revealed. The vertical pitch required for the 
project is 35.74 mm, centerline to centerline. The nylon 
bundle plates, as manufactured, had a pitch that was 
considerably less. In an effort to correct the tube mis- 
alignment problem, a low-powered laser was proposed to align 
the bundle tubes. A new nylon endplate was fabricated for the 
exit endcap, at the correct pitch with the tube bundle reduced 
to four vertical tubes,due to constraints in the endcap 
openings. An aluminum mount for the laser was fabricated by 
machine shop personnel that fit snugly into the tube 
penetrations of the nylon block. Several test projections 
were conducted on a plastic template fitted to the entrance 
endcap to minimize misalignment possibilities. The plastic 
template was scribed and cut with particular attention paid to 
correct pitch and used to manufacture the entrance nylon 
endplate. Subsequent visual inspection during experimental 
runs revealed that misalignments and bowing had _ been 


corrected. 
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co Vapor Superheat Problems 


During inate tad experimental runs, observed 
discrepancies in the measured temperature increases across the 
condenser tube lengths when compared to predicted temperature 
increases, prompted consideration of the possibility of vapor 
superheat occurring as vapor (at saturation temperature less 
than room temperature) flowed from the boiling chamber through 
the riser and up the sides of the condensing chamber. A two- 
pronged strategy was developed to minimize this possible 
effect. First, the apparatus shell was completely insulated 
with 12.7 mm thick foam insulation and all ancillary tygon 
tubing was insulated with double-wrapped, 3.18 mm thick foam 
insulation. Secondly, all experimental runs were performed at 
vapor saturation temperatures above the ambient temperature. 

4. Contamination Problems 

Visual inspection of the condensing tubes during 
operation, correlated with observed disparities in measured 
temperature increases across condenser tubes, revealed a 
structured, crystalline surface or matrix formation on the 
tube surfaces that inhibited heat transfer. The first 
appearance of the inhibiting matrix occurred during runs with 
smooth copper tubes and with the coolant inlet temperature at 
approximately -20 C. The working fluid, R-114, was 
transferred to the storage unit except for an oily residue 
that remained in the bottom of the boiling chamber. This 


residue was not analyzed, but was assumed to be machining oil 
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that remained after apparatus fabrication. The system was 
flushed twice, once with acetone and once with R-113, and the 
boiling chamber was scrubbed out. The smooth copper tubes 
were re-installed after a complete cleaning of their outside 
surfaces, and the experimental runs were repeated under the 
same conditions. Visual inspection again revealed the 
presence of the matrix formation on the tubes, and the 
possibility of water contamination in the R-114 was 
conjectured as a probable cause. The matrix was thought to be 
ice crystals solidified on the tube surface. The decision was 
made to continue experimental runs with a different 
refrigerant (R-113) and a different type of condensing tube 
(copper-nickel 1024 fpm low integral-fin tube), with the aim 
of understanding the conditions that caused the phenomena to 
exist. First, a run was made using R-114 with the copper- 
nickel finned tube under exactly the same conditions as 
described in the smooth copper tube runs. Visual inspection, 
during this run, revealed an apparent thickening of the fins 
at the top of the tube when compared visually with the second 
tube (see Figure 4.1(a) fin normal appearance and Figure 
4.1(b) fin's thickened appearance). This appearance of 
thickening gave credence to the belief that water 
contamination was causing an ice layer to form on the top tube 
between fins. The refrigerant was then changed to R-113 and 
the run was repeated, with the inlet cooling temperature at 


approximately -20 C. Again, a marked thickening of the finned 
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surface appeared on the top of the first tube. The next 
experimental run was conducted on the same finned tube with R- 
113, but at a coolant inlet temperature slightly above 0O C. 
In this situation, no visible thickening of the finned surface 
was observed, but disparities in temperature increases across 
tube lengths persisted. The run was repeated several days 
later to observe whether the phenomenon consistently repeated 
itself. Finally, the smooth copper tube was substituted for 
the finned tube and the experimental run repeated at the 
higher coolant inlet temperature. Noticeable surface 
contamination, having an orange peel appearance, recurred. In 
addition, during this run, there were indications that a 
pocket of non-condensable gases formed in the top of the 
condensing chamber. As every effort had been made to evacuate 
the apparatus prior to commencing the experimental run, and no 
evidence exists to support outside leakage of air into the 
system, the possibility exists that contamination within the 
system produced the non-condensable gases during boiling. 
Two possible sources of the contamination are 
conjectured. The first is that high temperatures on the 
heater tubes and poor circulation in the boiling chamber are 
combining to break down the refrigerant molecules producing 
some type of hydro-carbon. The second is that a chemical 
reaction is occurring between the refrigerants and gasket 


material in the apparatus producing a hydro-carbon. Time and 


Qe 


funding limitations have prevented further investigation of 
the problem. 
5. summary 

The original test plan called for the testing of 
various enhanced surfaces in a Simulated bundle during 
condensation of R-114. The contamination problem described 
above, coupled with the time delays inherent in achieving 
solutions to the other aforementioned encountered problems, 


severely limited the results of this thesis. 


B. DATA REDUCTION 
1. Description of Program Capabilities 

The computer program, DRP1F, that collects’) and 
processes raw data, in conjunction with a Hewlett-Packard 
computer/data acquisition system described in Chapter II, is 
listed in its entirety in Appendix C. The program is designed 
to calculate and plot heat transfer parameters for a variety 
of different tube bundles, utilizing R-114, R-113, or steam as 
the working fluid. The program has the added capability of 
allowing testing for single tube performance. The program 


consists of five main sections, as follows: 


1. Driver Progran, 
2. Main Progran, 
3. Property Subroutines, 


4. Modified Wilson Plot Subroutine. 


5. Plotting Subroutines. 
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Of these five sections, only the modified Wilson Plot 
subroutine will be described in lengthy detail. 

The driver program permits the user to take data, 
reprocess data, or plot re-processed data through various 
subroutines. The driver program is listed in lines 1000 
through 1125, in Appendix C. 

The main subroutine (lines 1130-2315) can be divided 
into five parts. The first part allows the user to select the 
physical parameters used in data reduction. The selection of 
parameters consists of the working fluid to be used (R-114, R- 
113, or steam), the vapor saturation temperature to be used 
(derived from averaged thermocouple readings in either the 
vapor section or liquid region of the apparatus), the 
instrumented test tube type, and whether data will be taken 
on a bundle or individual tubes. The second part allows the 
user to reprocess data with the same selection of physical 
parameters described above, but calls the modified Wilson Plot 
subroutine to calculate the inside and outside coefficients 
used in the correlations. Basic data reduction takes place in 
the third part (lines 1945-2125). Subroutines are called to 
calculate the properties of the ethylene glycol-water 
solution. These properties, with the physical tube 
dimensions, are used to calculate velocities and Reynolds 
numbers. Low coolant velocities and high viscosities prompted 
the use of twisted tape inserts (thickness 0.559 mm, with a 


pitch for a 180 degree twist of three times the tube's inner 
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diameter) to increase the inside heat flux. The inside 
Nusselt number with a twisted tape insert, provided by Hong 


and Bergles [Ref. 21], is given by the correlation: 
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where the Reynolds number for coolant flow was given by: 


Re_ = 4-T/nen+(D, - 46) (4.2) 


where, 6 is tape thickness. 
The outside heat transfer coefficient is then given by 


the well-Known summation of resistances to heat transfer: 
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which algebraically reduces to: 
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The heat transferred to the coolant (Q) is given by the 


relationship: 
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The heat flux (Q") is subsequently calculated by dividing by 


the outside surface area: 


pe! (4.6) 
Q A 
O 
where A, = ™D,°L. And, the overall heat transfer coefficient 


(U,) is given by: 


where the Log Mean Temperature Difference (LMTD) 1s defined 


to be: 
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It should be noted that the wall resistance due to fouling was 
assumed to be negligible for the purposes of calculation. 

The fourth part of the main subroutine creates a raw 
data file (lines 1680-1835), allowing subsequent reprocessing 
by the modified Wilson subroutine. The fifth part of the main 
subroutine provides a printed output both while taking initial 
data and subsequently after reprocessing data. 

The third major section of the computer program 
calculates fluid properties through called functions for both 


the working fluid and the coolant. The calculated coolant 
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properties as a function of temperature (lines 2330-2615), are 
kinematic viscosity, specific heat, density, Prandtl number, 
ema conductivity. 

The fifth major section of the program provides 
plotting routines for the output files generated in the 
program. The relationships graphically displayed are heat- 
transfer coefficient ratios (either based on the first tube in 
the bundle or as a ratio of the Nusselt value) plotted against 
tube position, heat transfer coefficient plotted against 
either the heat flux or temperature rise of the coolant, and 
the X-Y plot generated from the modified Wilson Plot results. 

2. Modified Wilson Plot 

The modified Wilson Plot, as outlined by Marto [Ref. 
22], accomplishes an indirect measurement of the outside heat- 
transfer coefficient. The implicit assumption is that the 
overall heat transfer coefficient (U,) 1s reliably known from 


data and, therefore a summation of heat transfer resistances 


is assumed. This summation relationship is given by: 
A 
he 1/2 4.9 
SES t RA t RAG tio (4.9) 
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where resistance due to wall fouling is assumed equal to zero. 
The summation equation 1s transformed to a linear 


relationship, as follows: 
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where: 
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which results in the simple linear form, Y = mX + b. ate 
should be noted that Theta (defined in line 3290) is derived 
from the Hong and Bergles [Ref. 20] relationship for the 


inside Nusselt number, given by: 
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Further, F (defined in line 3330), is derived from Nusselt's 
relationship for a horizontal cylinder subjected to a constant 


heat flux (Ref. 4], given by: 
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X and Y values are calculated from raw data and the data are 
fit with a least squares approximation. Initial assumed 
values are taken from the aforementioned correlations, and the 
solutions iterated to find the inside coefficient (C,) and the 


outside coefficient (C,) that fit the data, where: 





— ee ; = 
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the slope (m) is given by: 
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and the intercept (b) is given by: 
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The accuracy of this method relies heavily upon the number of 


data points taken and the range of velocities utilized. 
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V. RESULTS AND DISCUSSION 


A. EXPERIMENTAL SEQUENCE 


The data runs are summarized in Table 5.1 according to 


tube type, refrigerant used, and approximate operating 
parameters. Specifics for each run are provided in the 
following section under the run title. 
TABLE 5.1 

SEQUENTIAL LISTING OF DATA RUNS 
Run Tube Refrigerant Coolant Vapor 
Oe ll id y= Type Used T inlet T saturation 
SMT0O2 smooth R-114 a7. .C 1/7 .8 7 
CNFTO1 finned R-114 = 2 ene 13.36 
CNFTO2 finned jeg LAL = 2b € 64.4 C 
CNFTO3 finned Roe les -3 to +6 C 52.256 
CNFT0O4 finned Asad Gd =—3 CO +66 46.3 C 
SMTO3 smooth Rs =]. €O +676 56.7 oe 


B. EXPERIMENTAL RESULTS 


a 


copper 
temperature was 17.8 C. 


ethylene glycol coolant inlet temperature was 


while 


condensing R-114. 


Vapor 


Smooth Tubes with R-114 and Low Coolant Temperature 


SMTO2 was an experimental run made with four smooth 


Ambient temperature was 20 C. 
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saturation 


experiment was conducted following a cleaning of the boiling 
chamber, two flushings of the apparatus, once with acetone and 
ence with R-113, evacuation to 27.5 in Hg., and filling the 
system from an unopened cylinder of R-114. After filling the 
system, the apparatus pressure stabilized at 17 psig. The 
system pressure was lowered to 10 psig by opening the coolant 
flow through the auxiliary condenser. When the desired gage 
pressure was reached, the lowest flow rate was set through the 
instrumented tube bundle, and all three heater units were set 
for power levels corresponding to approximately 10 kW total. 
Gage pressure was maintained nearly constant by controlling 
coolant flow through the auxiliary condenser. The system 
required approximately 35 minutes from lightoff to reach a 
steady state condition, indicated by nearly constant gage 
pressure. The R-114 appeared clear at the commencement of the 
Sata un. By the time the apparatus reached a steady state 
condition, the contamination matrix was fully formed and 
visible on the top two tubes in the bundle. Data were taken 
at nine different coolant velocities, after the voltage 
indicator on the data acquisition system for the bottom tube 
thermopile appeared to reach a fixed value following each 
velocity change. Data for the tube bundle run (Figure 5.1) 
clearly demonstrates an observed uncertainty of about 7% in 
the overall heat transfer coefficient (U,) at velocities close 
to 0.9 m/s. This observed uncertainty, in contrast to the 


Maximum calculated uncertainty of 4.5% (see Appendix A), 
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SMT02--Tube Bundle Performance 
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indicates the effect of the contamination on the measurement 
of the overall heat transfer. The most significant result of 
the contamination is the addition of an un-quantifiable 
resistance to heat-transfer that has produced heat-transfer 
coefficients approximately half of those predicted by the 
Nusselt correlation. Further, it is apparent that during 
bundle operation, the contamination degrades the performance 
of the tubes in a graduated manner with the first tube being 
the most adversely affected and the third tube appearing the 
least affected. 

Upon completion of the bundle run, the flow in all of 
the instrumented tubes was shut off. Gage pressure was 
maintained at 10 psig by controlling coolant flow through the 
auxiliary condenser. The coolant inlet temperature remained 
approximately -20 C. The apparatus remained in this condition 
until the tube surfaces appeared to dry off, which took 
normally about ten minutes. Data were then taken at the same 
nine flow rates, as for the bundle, but with flow through only 
one instrumented tube at a time, starting with the top tube. 
Between velocity changes, approximately 10 minutes was allowed 
for the temperature changes in the tube to take effect. Data 
runs were made only after inundation from the previously 
tested tube was no longer present. The R-114 in the boiling 
chamber remained visually clear for all single tube runs. 
Again, as in the bundle data, data from the single tube runs 


(Figure 5.2) fall within an acceptable uncertainty band. 
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From the conditions under which these runs were made, it is 
noteworthy and logical that the effect of the contamination 
increases in magnitude as a run progresses. The first tube 
gives the lowest results possibly because of residual 
contamination from the bundle run. The first and second tubes 
produced approximately the same resultS as_ previously 
discussed in the bundle run, however the third and fourth 
tubes gave lower results when operating as single tubes than 
when operating in a bundle, possibly because condensate 
inundation that provides a rinsing effect in the bundle 
operation,is no longer present during single tube operation. 
2. Finned Tubes with R-114 and Low Coolant Temperature 

CNFTO1 was an experimental run made with copper-nickel 
low integral-fin tubes condensing R-114. Vapor saturation 
temperature was 18.3 C. The ethylene glycol coolant inlet 
temperature was at -21 C. The experiment commenced following 
the completion of the SMTO2 single tube runs, transfer of the 
R-114 to the auxiliary storage unit, installation of the 
copper-nickel finned tubes, evacuation of the apparatus to 
27.5 in Hg, and a refill of R-114 from the storage unit. This 
process took approximately six hours to complete. The same 
data taking procedure as outlined for the smooth tube runs was 
followed. The appearance of the R-114 was clear at the 
commencement of the run. By the time the system reached 
steady state, a marked thickening of the fin areas on the top 


tube when compared to the fins of the second tube was apparent 
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(see Figure 4.1). The R-114 appeared clear at the end of the 
bundle run. Data from the bundle run (Figure 5.3) show an 
increase in the overall heat transfer coefficient of 
approximately 100% when compared to the bundle results for the 
smooth tubes SMTO2 (Figure 5.1). The coolant velocity 
differences between the smooth tubes and the finned tubes is 
due to the smaller inside diameter of the finned tubes at the 
prescribed mass flow rate. There is no disparity between 
calculated and observed uncertainties given the calculated 
uncertainty band, but performance remains approximately half 
of that expected when compared to the values obtained with the 
Nusselt correlation and enhancement ratios reported by other 
investigations. Once again,the first and second tubes 
demonstrate the greatest degradation in performance due to the 
contamination. The third and fourth tubes either receive less 
contamination or derive enhancement from condensate inundation 
flow. From observations, the unknown contaminant had a higher 
surface tension than R-114, but to what extent the contaminant 
kept fin root areas flooded and subsequently negated the 
enhanced surface effect, remains obscured in the uncertainties 
of the measurements. 

The single tube runs followed the procedures outlined 
for the single tube runs of the smooth copper tubes. However 
at the end of each single tube run, a marked thickening of the 
fin surface at the top of the tube was noticeable. Vapor 


saturation temperature and coolant inlet temperature were the 
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CNFTO1L--Tube Bundle Performance 


Figure 5.3 


Same as in the bundle run. Data from the single tube runs 
(Figure 5.4) demonstrate an enhancement of approximately 76% 
in the overall heat transfer coefficient when compared to its 
counterpart in SMTO0O2. The spread in the data 1S approximately 
7% and, given the calculated uncertainty band, appears 
reasonable. Performance remains lower than expected from 
existing studies, and clearly the magnitude of the 
contamination effect increases as each run proceeded. 
3. Finned Tubes with R-113 and Low Coolant Temperature 

Upon completion of the single tube runs described in 
CNFTO1, R-114 was transferred to the storage unit. The 
apparatus was drained of residual refrigerant, the system was 
evacuated to 27.5 inHg, and filled with freshly distilled R- 
113 by using apparatus vacuum to promote flow from the R-113 
container. After filling, system vacuum was at 11 inches Hg. 
The system was brought to a gage pressure of 5 psig 
corresponding to a vapor saturation temperature of 64.4 C. 
Coolant inlet temperature remained the same as described in 
GNF POdt 7 ac. Vapor saturation temperature was monitored 
through the system pressure, and controlled by regulating 
coolant flow through the auxiliary condenser. The same nine 
Gata points described previously were taken in accordance with 
the stated procedure. The appearance of the R-113 at the 
beginning of the run was clear and clean. Upon the conclusion 
of the run, however, the R-113 had a slight yellowish tinge. 


During the run, the same marked thickening of the fin surface 
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CNFTO1--Individual Tube Performance 


Figure 5.4 


described in CNFT0O1, and shown schematically in Figure 4.1 
occurred on the top tube. The data from the bundle run 
(Figure 5.5) demonstrates only a slight improvement in the 
overall heat transfer coefficient when compared to the bundle 
Gata for SMTO2 (Figure 5.1), and in fact represents an 
approximate 80% decrease in the overall heat transfer 
coefficient when compared to the performance of the bundle in 
CNFTO1 (Figure 5.3). In view of the fact that these data were 
collected after transferring the R-114 to the storage unit and 
filling with freshly distilled R-113, it can only be assumed 
that the contamination was freshly created with the R-113 and, 
in fact is associated with the apparatus. It is also 
possible, based on these results, to infer the presence of 
non-condensable gases generated by the contamination, although 
no indications of non-condensable gases were detected during 
the run. In view of these facts, the reliability of the 
calculated values for the overall heat transfer coefficient is 
extremely doubtful even though data spread falls within the 
acceptable uncertainty band. 

Upon completion of the bundle run, data were taken for 
each of the tubes individually at the same system conditions 
described for the bundle run, and in accordance with the 
procedure described for single tube runs. The appearance of 
the R-113 at the completion of the single tube runs was 
slightly-darkened, with a yellow tinge, but otherwise appeared 


free from any floating contaminants. On each of the _ single 


50 


=. 00 






1.85 


LEGEND 


1.70 


\f> 
v) 


too 1.40 


VEL. (M/S) 


eG, 


TENE 


0.95 


0.80 


QO;O08 OUGA OODOLZ OOSGI VQOOY VOGSG QUOSG DAUGh YQoor AOdsSt VBYUE 
(eCard /M) 1 


0.65 


20 


ei 


CNFTO02--Tube Bundle Performance 


PieuEe. a. 


tube runs, a marked thickening of the fin surface developed in 
the course of the run. There exists reasonably good agreement 
between the single tube data (Figure 5.6) and the single tube 
Gata from CNFTO1l in the calculation of the overall heat 
transfer coefficient, and no disparity exists between the data 
spread and the calculated wincans sents band. The apparatus was 
evacuated for five minutes prior to beginning the single tube 
runs and if non-condensable gases were present during the 
bundle run, this explains the improved performance of the 
tubes during individual runs. 
4, Finned Tubes with R-113 and High Coolant Temperature 
This data run was made with the same copper-nickel 
tubes described in CNFT0O1 and CNFT02. The system was allowed 
to sit overnight at 11 in Hg, while the sump temperature was 
brought to an operating range of -3 C to +6 C. No change in 
the system vacuum over this 12 hour period occurred. The 
decision to run the system at a higher coolant temperature was 
motivated by a belief that the contamination of the system was 
water, and that at the colder coolant temperature, an ice 
sheath was forming on the tubes inhibiting heat transfer. The 
appearance of the R-113 before commencement of the bundle run 
was a medium yellow tinge. The system was brought to a steady 
state pressure corresponding to a vapor saturation temperature 
of 52.2 C. Difficulty was encountered in maintaining positive 
system pressure with the auxiliary condenser operating. The 


auxiliary condenser was therefore turned off and steady state 
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waS maintained by slight power adjustments to the heating 
elements in the boiling chamber. Nine sets of data were taken 
in accordance with the procedure for bundle runs previously 
described. No observable thickening of fin surfaces during 
the run occurred. The data from the bundle run (Figure 5.7) 
shows a slight increase in overall heat transfer coefficient 
over the bundle results obtained in CNFTO2 (Figure 5.5) and 
this is possibly explained by a decrease in the viscosity of 
the contaminant when exposed to the warmer tube surface. 
Single tube runs were made on each of the tubes, in 
accordance with procedures outlined previously, at the same 
coolant inlet temperature as in the bundle run, but at a 
system pressure corresponding ELS) a Vapor saturation 
temperature of 64.4 C. No change in the appearance of the R- 
113 was detected. Again no thickening of the fin surfaces was 
detected. The single tube results (Figure 5.8) demonstrate a 
corresponding range of values in the overall heat transfer 
coefficient when compared to the results for single tubes in 
CNFTO2 (Figure 5.6), but is significantly higher than the 
bundle results. It is possible that non-condensable gases 
were present in bundle operation, but were eliminated by 
evacuation prior to single tube data being taken. No 
indications of non-condensable gases were detected during 


bundle or single tube operation. 
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5. Finned Tubes with R-113 and High Coolant Temperature 


This data run was made in order to see if the results 
produced in CNFT03 were repeatable. The system operating 
pressure corresponded to a vapor saturation temperature of 
46.3 C, and the coolant inlet temperature remained unchanged 
from the coolant temperature in CNFT03. Nine data sets were 
taken in accordance with the bundle run procedures outlined 
previously. No thickening of fin surfaces during the run was 
observed. The color of the R-113, at the completion of the 
run, waS amber. The bundle data (Figure 5.9) does demonstrate 
that the results of CNFT03 were repeatable at least within the 
calculated uncertainty. 

Single tube runs were made in accordance with the 
procedures discussed previously. Nine data sets were taken on 
each tube, at the same conditions specified in the bundle run. 
No thickening of fin surfaces was observed during the runs, 
and the R-113 color remained amber. There is no major 
disparity between the single tube results (Figure 5.10), and 
the single tube results for CNFT03, given the calculated 
uncertainty band. The single tube results were significantly 
higher than bundle performance, and the build up of non- 
condensable gases during bundle operation and evacuation of 
the apparatus prior to commencing single tube runs is 
conjectured as a possible explanation. At no time during 
bundle or single tube operation were indications of non- 


condensable gases detected. 
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6. Smooth Tubes with R-113 and High Coolant Temperature 


SMTO3 was a data run made on the same aforementioned 
copper tubes, but with R-113. Vapor saturation temperature 
was 56.7 C. Coolant inlet temperature ranged between -1 C and 
+6 C. The experiment was conducted following a test run with 
R-113 and finned tubes and subsequent replacement of the 
finned tubes, and evacuation to 11 in Hg. While the tube 
change was being made, the R-113 remained in the boiling 
chamber exposed to the atmosphere for a two hour period. The 
R-113 appeared light amber in color at the commencement of the 
bundle run. The desired gage pressure wasS maintained by 
controlling the boiling chamber heaters through their 
corresponding variac controllers. The auxiliary condenser was 
not utilized because the additional coolant flow provided by 
the auxiliary condenser made it difficult to maintain the 
desired system pressure. During the data run, the top two 
instrumented tubes assumed an orange peel textured or dimpled 
appearance. The R-113 had a darker amber appearance at the 
completion of the bundle run. Again, data were taken at the 
same nine coolant flow settings following the procedure 
outlined previously. Data taken during bundle operation 
(Figure 5.11), while falling within the acceptable uncertainty 
bands, clearly demonstrates the presence of non-condensable 
gases as well the presence of the contamination acting to 
inhibit the overall heat-transfer performance. The top tube's 


very poor performance can be attributable to its location in 
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a pocket of non-condensable gases. And the bundle performance 
as a whole, when taken in comparison with the performance 
demonstrated in the bundle run of SMT02, indicates an 
increased concentration of the contamination and the presence 
of non-condensable gases. 

Upon completion of the bundle run, coolant flow 
through the instrumented tubes was shut off and gage pressure 
maintained by controlling flow through the auxiliary 
condenser. While no coolant flowed through the instrumented 
tubes, the textured appearance noted during the bundle 
operation gradually transformed to small droplets on the top 
two tubes which then coalesced into larger ellipsoid shaped 
droplets that were held on the tube by surface tension forces. 
Data was taken on the first two tubes in the bundle following 
the procedure outlined previously for single tube runs. At 
the completion of the second run, it was noticed that the top 
view port flanges were cold to the touch while the bottom view 
port flanges were warm to the touch. The assumption was 
therefore made that non-condensable gases had collected in the 
top of the condensing chamber. Before evacuation of the 
apparatus could be effected, the pyrex glass in the top center 
view port cracked. As a result, the apparatus was secured and 
no further data taken. It should be noted that, prior to the 
view port breaking, there was no evidence to suggest that the 
apparatus was other than air tight, so that the source of non- 


condensable gases was most probably due to their generation 


G2 


within the boiling chamber. Data taken in the single tube 
runs (Figure 5.12) reflects not only the effects of the 
contamination, but also the effects of the presence of non- 
condensable gases as well. This data showed the greatest 
variance between calculated and observed uncertainties in the 
overall heat transfer coefficient. This data is limited to 
the performance of only the top two tubes, due to the afore- 
mentioned breakage of the view port and the presence of non- 


condensable gases. 


C. SUMMARY OF RESULTS 

The comparison between finned and smooth tubes shows an 
enhancement ratio in the overall heat transfer coefficient for 
the finned tubes compared to the smooth tubes of approximately 
ere © While this enhancement in the overall heat-transfer 
coefficient might translate to a vapor-side heat-transfer 
coefficient enhancement ratio of approximately 4.0 to 4.5, 
based upon the calculated values for the vapor-side heat- 
transfer coefficient for the tested smooth tubes when compared 
against the Nusselt's correlation for smooth horizontal tubes 
and the reported enhancement ratios from other finned-tube 
investigations [{Refs. 5-13}, it 1s estimated that the effect 
of the contamination during these tests has been a degradation 
of up to 50% in the vapor-side heat-transfer coefficient. The 
contamination has resulted in the promotion of another 
resistance to heat transfer that has no predictable 


characteristics and can not be quantified. This contamination 
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must therefore be removed before any quantifiable heat- 


transfer results can be obtained. 
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A. 


VI. CONCLUSIONS AND RECOMMENDATIONS 


CONCLUSIONS 


Based on the data gathered during this investigation for 


condensation of refrigerants R-114 and R-113 in the multi-tube 


apparatus, the following conclusions are reached: 


es 


B. 


The construction of a multi-tube condensation test 
apparatus begun by Murphy [Ref. 3] and Zebrowski [Ref. 
4) has been completed, instrumented and integrated with 
the various support systems. 


Data reduction software has been developed that produces 
high confidence results from experimental readings. 


Data were collected from condensation experiments with 
R-114 and R-113 on smooth and finned tubes demonstrating 
an enhancement ratio of approximately 2.0 in the overall 
heat transfer coefficient for the finned tubes over the 
smooth tubes. 


Serious apparatus internal contamination has degradated 
the value of the overall heat-transfer beyond acceptable 
limits, making further data reduction meaningless. 


RECOMMENDATIONS 


Based upon the data gathered during the operation of the 


multi-tube condensation test apparatus, and the experience 


gained in construction of thiS apparatus, the following 


recommendations are made: 


1. 


Samples should be taken from the R-114 and R-113 used 
in the experimental runs and a chemical analysis 
performed on these samples to assist in determining the 
source of the contamination in the system. 


Modifications to the apparatus, as dictated by the 


results of the chemical analysis of the refrigerants, 
should be made. 
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New flow meters allowing a greater range in coolant 
velocities should be installed. The greater range in 
coolant velocities will enhance the accuracy of the 
modified Wilson plot method in calculating the outside 
heat-transfer coefficient. 


The temperature control system for the eight ton 
refrigeration unit should be re-evaluated to determine 
if tighter control of sump temperature can be 
maintained. 


After appropriate modifications to the apparatus have 
been made, the smooth tube and finned tube experiments 
should be repeated and contrasted with the results 
presented here, in order that baseline tube performance 
might be established. 


During operation of the system, flow should be throttled 
to the coolant supply pumps with the valves on the 
coolant inlet side. This will increase the pump life 
and reduce noise in the laboratory space. 
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APPENDIX A 


UNCERTAINTY ANALYSIS 


A certain amount of uncertainty exists in any engineering 
measurement. These uncertainties arise both from known 
sources, such as calibration and measurement errors from 
sensing devices, and also unknown sources, such as operator 
experience and unexpected experimental reactions. The 
procedure used in this uncertainty analysis is based on the 
Kline-McClintock [Ref. 21] method. This method assumes that 
a result R is a function of the See et SL@e that contributewes 
that result. Therefore the uncertainty in a result is a 
function of the uncertainties in each of the variables. This 


method is expressed mathematically by: 


oR auch ae ee ORS pe 262 172 
RT SE >? 1 ‘V; Rn age oe (A.1) 


The measurement of the overall heat transfer coefficient 


(U,) performed in this thesis is given by: 
U = Q (A.2) 


oO A -LMID 


where the heat transferred (Q) is given by: 


Q = TC - (AT) i 
P 
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and the Log Mean Temperature Difference (LMTD) is defined as: 


Therefore the uncertainty in the measurement of the overall 


heat transfer coefficient (U,) is given by: 





: 6C 6D 
2 z or, 2 Zz Guene, 2. 1/2 
se = 1? + Pe A s 2 AY (A.5) 
O : p O 
i 
6A 6D 
where —~ reduces to = , and the uncertainty in the LMTD is 
O O 
given by: 


Ltn 7.2 2 BN GD 


where A is given by the expression: 


ZS OT 2 De rere, 
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where B is given by the expression: 
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where C is given by the expression: 


69 


6Tc i 
@=  — 2. — aaa 
(Tat TC) 1 satis dy 
O” lip -Tc 
sat oO 


It should be noted that the uncertainties for mass flow 
rate (Pr) and the coolant temperature rise across the 
condensing tubes (AT) were calculated during calibration of 
the flow meters and thermopiles respectfully. Reference 
material providing data on thermophysical properties, 
specifically the specific heat (C,), listed no uncertainties 
associated with the curves for aqueous-ethylene glycol 
solutions. The composition of the aqueous-ethylene glycol 
solution by weight percentages was calculated by measurement 
of the specific gravity and the change in this composition 
over the experimental period was determined to be negligable 
and therfore this uncertainty has been ignored. 

The uncertainty in the overall heat transfer coefficient 
represents the maximum uncertainty present from the 
calibration of the four flow meters and the maximum 
uncertainty in the calibration of the four thermopiles, 
regardless of association, and therefore representS a 
conservative estimate of the uncertainty present in each 
calculation of the overall heat transfer coefficient (U,). 
This maximum uncertainty was calculated to be 4.5 percent. 

The disparity between the calculated uncertainties and the 
uncertainties observed by the spread in data particularly in 


the single tube runs, can be accounted for by the as-yet 


7.8) 


undetermined chemical or phase reactions generated by the 
contamination in the system. 

The values of the uncertainties in the measured variables 
that were used to calculate the uncertainty by the Kline- 


McClintock method are listed in Table A.1. 


TABLE A.1 
UNCERTAINTY VARIABLES 


VARIABLE VALUE REMARKS 


FE 
iy 
O 
a) 


From flow meter 
calibration data 


6C 
==" 0.00 Not available 
p 
7 Oe.0 ii From thermopile 
calibration data 
6D 
_ o° OOo From micrometer 
Dy name plate data 
OLMTD Listed below by 
LMTD run title 
SMTO2 62-005 
CNFTO1 0.002 
CNFTO2 0.001 
CNFTO3 0.003 
CNFT0O4 O2,0:02 
SMT0O3 O20 7 


yl 


APPENDIX B 


LIGHTOFF AND SECURING PROCEDURES 


fave SYSTEM LIGHTOFF 

1. Push the starter button in the control box for the re- 
circulation pump. This control box is located on the bulkhead 
above the re-circulation pump in the outside area adjacent to 
the refrigeration unit. 

2. Turn the switch on the refrigeration unit control 
panel, located in front of the refrigeration unit to the 
"auto" position after passing through "on" position. 

3. Set the desired temperature on the roughly graduated 
Fahrenheit scale on the control panel thermostat. It requires 
approximately four hours to chill the sump 40 degrees C. The 
thermometer located on the side of the ethylene glycol sump 
must be monitored to ensure the desired sump temperature is 
attained and maintained. Slight adjustments in the 
refrigeration unit thermostat can be expected due to the 
coarseness of its scale. 

4. Energize the heater variacs by switching on the 
breakers in the breaker panel located in the laboratory space 
on the bulkhead next to the counter. 

5. Set the heater variacs to the desired position, after 
ensuring that the switch panel for the heater tubes located on 


the apparatus has all switches in the "on" position. Monitor 


Ue 





apparatus pressure through the pressure gage at the top of the 
apparatus, ensuring system pressure does not exceed 30 psig. 

6. Turn on the pump motors by pushing down on the arm of 
the appropriate breaker box for the pumps located on the 
bulkhead next to the ethylene glycol sump. The pumps are 
marked "auxiliary condenser" and "instrumented tube 
condenser," respectively. Flow in the auxiliary condensate 
system can be controlled with the individual gate valves 
located at the coil penetrations on the apparatus. The 
auxiliary condenser will produce the fastest adjustments to 
system pressure if pressure is rising too quickly. Flow 
through the instrumented tubes can be controlled by the ball 
valves located at the bottom of the respective flow meter. 

7. Throttle down coolant flow to the supply pumps with 
the valve located before the pump suction. This will increase 


the pump life and reduce noise in the laboratory space. 


B. SECURING PROCEDURES 

1. Turn all variacs to the zero position and switch off 
all breakers in the power panel on the bulkhead. 

2. Turn the breakers for the pumps to the off position at 
the switch boxes near the ethylene glycol sump. 

3. If apparatus will not be operated for an extended 
period, turn the switch on the refrigeration control panel to 
the "off" position after passing through "on." 

4. Allow the re-circulation pump to operate for at least 


ten minutes after switching off the refrigeration unit to 
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dissipate any back pressure in the system; then secure the 


pump. 
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APPENDIX C 


DATA REDUCTION PROGRAM 


DRPIF 
SE: This program collects and processes condensation data for 


the R-114 tube-bundle apparatus. 
ED- NOVEMBER 2, loos 
ED: November 8, 1988 


CHANGE HOD=3 AND USING CALCULATED CI FROM WILSON 


ant jes 1 

WSNe 4x  SGLECH OPTION!” 

USING "6x,""@ TAKING DATA OR REPROCESING PREVIOUS DATA”*” 
USING “6X,""1 PLOTTING H VS DELTA-T""" 


Voom « PLOTTING HRAgeUS N° >? 
USING "bx 3S PLOTTING WILSON” 
USING 6x 4 uPURGE Filles: ’ 

DSING “6%, ° 5S XYREAD ”* 

PSING Ox, 96 NUSSELT ESTIMATE 
ER 1S 701 

Icall 


all=@ THEN CALL Main 
all=1 THEN CALL Plot2 
all=Z THEN CALL Protl 
el l=27 THEN CALL Plot3 
all=4 THEN CALL Purge 
all=5 THEN CALL Xyread 
all=6 THEN CALL Nusselt 


ain 
Bes 307) 
Pld? Sl ft 


COM /Nus/ Iin,Tsat.Qdol.Hnus 


com (7 
Eom 7 
BBO ures 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
READ 
DATA 
READ 
DATA 
DATA 
READ 
DATA 


Will/ Doa(4),D18a(4),Kma(4),lact 

Wil2/ Delta,Isat,Nsets,Hod,Cia(3),Alpaa(3) 
memes). Fm aaeemt( $8). 10'5).1C 8), Hot 5 G00 5), U0l 35) 
Goon 2Ovod. 4,34.5223- 
Ado. 25,8 1 woo, 26: 
Mol ole 46. 4.46.6 .54. 
S770. 55,505.5,55.6,58. 
Baeewe0 boca yo7., 7,445 
U2. 7 05,6562, 69.7548. 
80.6,50,76.3,76.7,55. 
88.4,55,83.3,83.8,60. 
96.3,60,90.2,90.8,6S. 
Fma(e) 

Bee? 2 ,S5.172,5.172,5.172 

Cia(*) 
VelO@cbUSie75727.94369,-75/345.6295, 76025595...81 
=9247486589,.6.97688E11,-2.66192E13,3.94078E14 
Cle) 

@.015875,20.014000,2.20,0.0,0.2 


OmNnNnNan rk wu 


ihe 


VATH 0.019259,0.0101b0,0.0,0.0,U.u 
DATA 386.0,42.975,.0.0,0.0,0.0 
DATA @.0005588, 3 


READ Doa(e),Dial(*#),Kmat*),Delta,Hoc |'| Hod=H/D: 
L=1.2192 ! Condensing length 

Jset=0 

BEER 


INPUT “ENTER MONTH, DATE AND TIME (MM:DD.HR-MM:S5)" Digs 
OUTPUT 709% 10> ete: 
Been 
INPUT "SELECT OPTION (@=DAQ, 1=FILE)",Im 
Ihard=1 
BEEP 
INPUT "WANT A HARDCOPY PRINTOUT (1=DEF=YES,Q0=N)",Ihard 
BEEP 
INPUT “SELECT (Q=R-1 14 )-STEAM. 2-R-11535-bo. 3 at 
Iin=1 
Isat=2 
BEEF 
INPUT “SELECT SAT TEMP MODE (O=LIG fp =UR oe Ele eer 2 Hee ee 
IF IThard=! THEN PRINTER IS 701 
IF Im=@ THEN 
Peas 
INPUT "GIVE A NAME FOR THE NEW DATA FILE" ,Files 
CREATE BDAT File$,2Q 
ASSIGN @6File TO Files 
BEEP 
INPUT “ENTER TUBE CODE” ,Itube 
Beer 
INPUT "ENTER EG CONCENTRATION (WT PERCENT)" ,Egrat 
ENTER 7@9:D0tgS 
OUTPUT ®File:Dtgs 
OUTPUT @File:ltube,Eorat, 0d! ,Od2 Des tes vEes 
Tact=1@ 
BEEP 
INPUT “SESECT <0=(CP 1=SEcONtne ==: , 1@=BUNDLE=DEFAULT >)" ,Iact 
PRINT 
PRINT USING “J0xX.”° “FILE NAME: VP Alea S 
PRINT 
FLSE 
BEEP 
INPUT “ENTER NAME OF EXISTING FILE" ,Files 
BEEr 
INPUT “ENTER NUMBER OF DATA SETS STORED" ,Nsets 
Iwil=] 
BEEF 
INPUT "WANT TO CALL WILSON? (1=DEFAULT=YES,@=NO)" ,Iwil 
IF Iwil=1 THEN 
BEEP 
INPUT “WHICH TUBE (@=TOP,1=SECOND,..,1Q=BUNDLE)",lact 
CALL Wilson 
IF Iharcg=1 THEN PRINTER 1S 721 


PRINT 

PRINT USING “10X,""“FILE NAME: mm”  i2h” Fakes 

PRINT 

IF Lact oe 
PRINT USING “10X,"~ "INSIDE COEFFICIENTS: “"s4\00.30-°2%5 | Craaee 
PRINT USING "1@X,""ALPHAS: Mea Ik, MZ. SOL. 2k) eee 
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@ai*) 


1S35 Bese 

1540 Pein) UstNowler. -INSIOE COEFFICIENT FOR TUBE *",0, ” 
Meeetact ,Cial lect) 

1545 PRINT USING “10X,°"ALPHA FOR TUBE “",0,°": 
D";Ilact Alpaa(lact) 

1552 ENO IF 

i eysis] EES 

is6s INPUT “ENTER CI VALUES (OEFAULT=S.172)" ,Cia(e) 
1572 END IF 

S75 ASSIGN @File TO Files 

1582 ENTER @File;Otgs 

Ue ss ENTER @File:Itube,Egrat,Odi,0d2,0d2,0d4,Dd5 


1S86 INPUT “ENTER TUBE COOE",Itube 

1592 END IF 

1595 Iout=l 

1602 BEEP 

16@5 INPUT "WANT TO CREATE AN OUTPUT FILE? (1=DEF=YES,Q@=NQ)",Iout 
1612 IF Iout=1 THEN 


ius its BEEP 

1622 iheOT SS ENTER Se NAME POR OUTPUT FILE ,Fouts 
1625 Cree 60Al Fouts,S 

163@ PoclGh CFout 10 Four? 

[64- BNO IF 

1642 |! 


1645 Oo=Doa(Itube) 

1550 01=O1a(Itube) 

1655 Km=Kma( Itube) 

1662 Axe=Pl*O0:°2/4 ' Cross-sectional area 
1665S Ao=PI*0o0+#L 

1672 Rm=Do#L06(00/D2 )/(2#Km) 


tees: | 
1682 IF Im=@ THEN 
168S PRINTER IS 1 
1696 BEEP 
1695 PRINT “SET FLOWMETER REAOINGS CORRESPONOING TO: “ 
1722 PRINT * ";Fma(Jset,1):"% OF METER 2 ANO HIT CONTINUE" 
1705 PAUSE 
1710 PRINTER IS 701 
1715 OUTPUT 709;"AR AFQ@ ALE VRS" 
1722 Nend=8 | INCREASE TO 9 IF FIVE TUBES IN BUNOLE 
eS FOR I=@ TO Nend 
1730 OUTPUT 709: "AS SA" 
IEG 21S) Vsum=0 
1742 FOR J=1 TO § 
1745 ENTER 7@09;E 
1752 Vsum=VsumtE 
i7Ss NEXT J 
1762 Emf( I )=#Usum/S 
1765 NEXT I 
1772 OUTPUT 709: "AR AFZ2Q@ ALZ23 VRS" 
77S FOR I=@ TO 3 
1782 OUTPUT 76S; "AS SA” 
1785 Vsum=2 
1792 FOR J=1 TO SQ 
1795 ENTER 7@9;E 
1822 Vsum=Vsumt+E 
18@S |! WAIT .25 
_ 1810 NEXT J 


1 


OD 3 


mO0.3 


1815 Tpo( I )=Vsum/S@ 


1822 NEXT I 

1825 EIESE 

1832 ENTER @File.Fm(*#), Emf(*+),Tp(?) 
1835 ENO IF 

184Q@ |! 

1845 '§ DATA ANALYSIS 

18S5@ |! 


18SS Nend=8 

1862 FOR I=@ TO Nend 

1865 TCI )©FNTvsv(Emf(l)) 
1872 NEXT I 

reds Ivease (1 (0410 ltt 297 
1880 ibrg=( 10394114) 072 

1892 Tsat=Tlig 


1835 END IF 
1896 IF Isatel THEN 





1902 Tsat=elvap 

19905 ENO IF 

19@6 IF Isat=2 THEN 

193Q7 Tsat=(Tvapt+Tlig)/2 


19@€& END IF | 
191@ FOR I=@ TO 4 


1 shUS Fm(] )=Fma(Jset,I) 

192@ NEXT ] 

ies JseteJset+] 

1930 PRINT 

1935 PRINT USING “1@X,""Data set number = ““,OD"; Jset 


1942 PRINT 
1945 Ibeg=2 


195@ lend=3 

1955 [Fo tact 10 1REN 
196@ Ibeg=lact 
1396S lenc=lact 


197@ ENG oe 
1345 FOR I=Ibeg TO Iend 


1380 Grad=FNGrad(Emfi1+S)) 

P3e5 Delt=ABS(Tpi1)/(Grad+le@)) 
1992 Tavg=I(1+S)+Delt+.5 

i345 Rhnoeg=FNRhoeg( Tavg,Egrat ) 
LU2e Nueg=FNNueg( Tavg,Eorat ) 

c@2S Mueg=Nueg*Rhoeg 

2@1@ Cpeg=FNCpeg( Tavg,Egrat ) 

eOrs Keg=FNKeg( Tavg,Egrat ? 

2022 Preg=Cpeg*Mueg/kKeg 

2025 Mdot=FNFmceal(1,T(I+S),Fmil)) 
2030 Veg=Mdot/(Rhoegs#A». ) 

Males Reeg=Veg*Di /Nueg 

~@42 Res=4eMdot/(Pl*Mueg+(D1-4*#Delta) ) 
2045 Qdot=Mdot+Cpeg*Delt 

2250 Qdp( I )=Qdot/fio 

2005S IF I=@ OR I=lIact THEN 

2262 Qdpl=Qdp(1) 

~@6S CALL Nusselt 

2872 END IF 

2@75 Latd=Oelt/LOG((]s@t—-1(1+S))/(iset—i(1+S)-werrs 
2082 Vo(l)=Qdp(I)/Lmtd 

2@6S IF Reeg< 4002 THEN 
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“2290 Nueg=Ciail )#(1+5.484E-3*Preg*.7*(Res/Hod)*1.25)".5 

2@95 E\LSe 

2180 BEEP 

2105 PRINT USING "1QX," "INCORRECT TURBULENT CORRELATION""" 

2118 Nueg=.@27+#Reeg*.8*Preg’ .5333*Cfeg 

2115 END IF 

212 Hi=Nueg*keg/D1 

cizs Ho( I )2#1/¢(1/Uo1 I )-Do/(DieH: )-Rm) 

2 WO IF I=Q@ OR J=lIact THEN 

2135 PRINT USING "10X,""Mass flow rate eo M727 oOe «Nedet 

2 ib PRINT USING "10X,"“Inside Tube Dia. (m.) eee MN OOE «bare. | tube 
) 

217 PRINT USING “10X,""18@ DEG OVER Dia. (HOD) = °°,MZ.3DE"; Hod 
2140 PRINT USING "1@X,""Inlet temperature ateemri? . SDE. ant 1+5 >) 
2145 PRINT USING "10X,""“Saturation temp (Deg C) = “",MZ.3DE";Tsat 
2146 PRINT USING "10X,""“DELTA Tape Thickness = ee ODE} OST ta 
2147 PRINT USING "10X,""DELT temp Oif. = = ieee” 7 Dett 
2148 PRINT USING "10X,°"Log. Mean temp O1f. Sac oor 7 lmto 
2150 PRINT USING “1@X," "Heat flux = eo weOe oops 1) 
2191 BoiNy USING ‘IUx, Conductivity E.6. zee” M2.00E i Keg 
Zoe PRINT USING "10X,""Conductivity Tube Metal = “",MZ.3DE"iKm 
2155 PRINT USING “10X,""°Prandtl number = “" MZ.3DE";Preg 
216@ PRINT USING "10X,""“Reynolds number Se SOEs neeG 
Age PRINT USING “"10X,""Reynolds number H&B S " MZ.3DE";Res 
2TES Pent USING ==10x%, Inside htc = "" M7Z.3DE";H1 
2166 PRINT USING “10X,""Inside NUSSULT NO. = ee oO bee NOeO 
2167 PRINT USING "10x," “OVERALL H.t.c. (Uo) Speen oe De” 7 Uo | ) 
2170 PRINT 

2S PRINT USING “10x, “"TUBE FM VEG DELT Uo 

HO HNUS""" 

2182 PRINT USING "1@X,"" # CL) (m/s) (K) Wm os 
ko 

Z Nels PRINT USING. “10x 20.4%, 20-08. 3% 72.00 .4% bie sie ie SDE 3X MZ. SDE 
eet ODE, SX 1 1+1,Fmél),Veg,Delt,Uo(l),Hofl),Hnus 

2192 SL Sle 

2s PRe MUS UNG * 10.530, 4x% , 30,00, 3%, 2.0004. OU. S02 tee OE BY MT? SOE 
Mee t+ 1 Fmt 1) Veg,Delt,Vot(l),Holl) 

2220 END IF 

2205 NEXT I 

mae) IF Im=@ THEN 

are Okaccpt=1 

2aee BEEF 

2ocS INPUT “OF TO ACCEPT THIS SET (1=DEFAULT=YES, @=NO)?" ,Okaccot 

Eee Meeomeccpt =. IHEN OUTPUT OFile; emt +) Emit*). tox ) 

cas END IF 

2-40 IF (Okaccpt=1 OR Im=1) AND Iout=l1 THEN 

2245 FOR I=@ TO 2 

2250 OUTPUT @Fout.Ho(!I),Qdpil) 

2255 NEXT I 

2260 END IF 

2265 IF Im=Q THEN 

a5 (4) Okrpt=1 

Len > BEEP 

2280 INPUT “WILL THERE BE ANOTHER RUN (1=YES=DEFAULT,Q=NO)",Okrpt 

2ee85 PevOkret=1 THEN 168@ 

2292 jal SIs 

2ess IF Jset<Nsets THEN 1680 

2202 END IF 

2305 ASSIGN @File TO * 


a2 


1210 
eae 
Zoa0 
25a5 
2550 
2395 
2348 
2345 
Ze50 
2250 
2360 
Za55 
2370 
Zone 
2380 
25o5 
2530 
Zo35 
2480 
2485 
2418 
2415 
2428 
2425 
2438 
2255 
2442 
c44s 
e452 
2455S 
2462 
2465S 
2472 
2475 
2482 
2485 
2490 
2455 
2500 
2505 
2510 
rey es 
2520 
Z525 
2550 
Z509 
2540 
2545 
Zoot 
2550 
2560 
2509 
Past) 
2505 
258@ 
Zoco 
2530 
2509 
2688 
facials 


IF Iout=1 THEN ASSIGN @File TO * 

SUBEND 

DEF FNGrad(T) 

Grad=-3.8778S7E-S-2*4.71428S7E-8eT 

RETURN Grad 

FNEND 

DEF FNKcu(T ) 

OFHC COPPER 252 TO 302 K 

TREeT+e75.15 

mh She NZS 

RETURN K 

FNEND 

DEF FNNueg(Tc,Egr) 

RANGE OF VALIDITY: -20 TO 2@ DEG C 

Tk#Tc+273.1S 
Nut=7.1196SQ7E-3-Tke(7.4863347E-S-Tk*#(2.6294943E-7-1k #3 .O83535227E-iae 
Nu2#4.9237638E-3-Tke( 4.921391 2E-S-Tke(1.6437S34E-7-Tkei .8333331E-10)) 
Nu3=8 .6586293E-3-Tke(8.8837902E-S-Tke#( 3,0495032E-7-Tk#3.4999996E-10)) 
AZ=(Nu3Z-2*NuZ2tNul )/280 

Al=(Nu2-Nul-940#AZ2 )/10 

AQ@=Nu1-42*A1-1764#A2 

Nu=A@+tEgr#(AlL+Egre#AZ ) 

RETURN Nu 

FNEND 

DEF FNCpeg(Tc,Egr) 

RANGE OF VALIDITY: @ TO 2@ DEG C 

fetes 2¢o. 15 

Upt=!.5 (VISSUE+STIkK eo. 3 

Cozi 4748 l25e+S+ih e025 

Cp3=9.S8G@SOGE+2+Tk#7.3 

A2=0 Co5-2*Cpzt+Cpl 7280 

Al=(Cp2-Cp1-988*AZ2 )/12 

A@=Cp1-40+*A1-1600¢Az2 

Cp=A@+Egre( Al+Egr Az ) 

RETURN Cp 

FNEND 

DEF FNRhoeg(T.Egr ) 

Ro lel GS0/7@S3E+5-14* 0s. 7G51285E—-14 1 €4..0637 le 3b 5) 
Rol=1.0748°72E+3-T*#(4.4266!19SE-1+T+4.0939706E-3 ) 
Ros=™!] .G885934E+3—-14¢5 .7S55655E-!+ Tea l2ersoSe—2) 
AZ=(Ros—-2*RoztRoel 1/700 

Al=(Ro2-Rol-98@8#Az 1/18 

AZ=Ro1-4@0*A1-160@*A2 

Ro=AQ+Egre(AlL+EgreAz ) 

RETURN Ro 

FNEND 

DEF PNEeee (Tt Eee) 
Pr=FNCpeg(T,Egr )*FNNueg(T,Esgr )*FNRoeg(T,Egr )/FNKegiT,Egr ) 
RETURN Pre 

FNEND 

DEF FNKeg(Tc,Egr ) 

RANGE OF VALIDITY: -2@ T0 28 DEG C 

TR=T¢G+2 755 G 
Ki=e2.2824708E-14+Tk#(5.5989286E-4+Th*#3.S5714°86E-7 ) 
K2=2,. S646616GE-1+T eC 25 Sod7eS 7 lE-dtihe 7s leees7 le 7 
K3=3.2138932E-1-Tke(3.00428S7E-4-Tke1.428S714E-6 ) 
A2Z=(K3-2*KZ+K1 1/280 

Al=(K2-Ki-900#AZ )/10 

AQ=K1-4@*¢A1-1600#AZ 

K#AQ+Eor*(Al+Egr Ac ) 


&0 


2612 RETURN k 

2615 FNEND 

PScec SEF FNTann(s ) 

file Fils P=EXP(» ) 

26352 Q=1/P 

2635 Tanh=(P-0@)/(P+Q) 

7642 RETURN Tann 

2645 FNEND 

2650 DEF FNTvsviV) 

2655 COMeAzCe, C7) 

2662 T=0(Q) 

2665 FOR I=) 10 7 

2670 T=ieet rev] 

2675 NEXT I 

2682 RETURN T 

2685 FNEND 

2692 DEF FNBeta(T) 

2695 Rop=FNRho(T+.1) 

2700 Rom=FNRho(T-.1 ) 

2705 Beta=-c/(Rop+Rom)*(Rop-Rom)/.2 
2710 RETURN Beta 

2715 FNEND 

2720 DEF FNPsat(Tc) 

2725 !' 0 TO 8@ deg F CURVE FIT OF Psat 
2732 Nf=1 568% let 32 

2735 Pa=5.945525+Tfe( .15352082+Tfe(1.4840963E-3+Tf*S.6150671E-6 )) 
2740 Po=Pa-14.7 


2745S IF Pg>®@ THEN | +£P516,-"1n HO 
2752 Psat=Po 

Z7o5 EESE 

2760 Psat=P9#29.92/14.7 


2765 END IF 
Yara) RETURN Psat 


emt>D rINcINy 


2 *BO DEF FNFmcal(I,1,fm) 

27185 DI eee 4), 6S(4),M104),M2(4),M3(4 > 

2790 Gimme ce Soe o 1 /49E-5,-5 8105E-5,-5.5079E-4.4.125E-4 

2735 DOr eae OZ O/E-5,—-2.7S596E-5,-4.309135E-2,-3.109Z27E-5 
2096 Opigpmienoteemenmele sce. ,)]. 1 14E-2,-2,.86717E-3,2.0 

2BOO Dance e soem (OS206—5,1.95749E-5.1,92237E-3,2.S6E46E-3 
2305 OO ieee Sbcies, 2. Co@Moas, 2. Udalte—-5,2 .2453E-3,2.92127E-35 
2B@6 DOHieibIlemon cs (VO2ZEoo2,0012Z1E-35,2.2040E-3,0:.2 

2612 PEGs bdtee, Bar" ),65( #7 M1Ce) M204) M3) 


aot} LES <SeerneN 

ZEAS ZVeB IX TO4MICI erm 
Zeee Z2£B2(1)+M2(1)¢Fm 
2B21 oa i OY op) 

ZUG. C=Z1+S+#20 

2824 EESE 

Zez> ivepeCl +20 erm 
2B26 22-6541 )+MS(1)*Fm 
|S Ra S=ee27 17-16 

2828 C=Z1+S¢7 


2830 ENG ie 

28635 Mdot=C+SeT 
2840 RETURN Madot 
2845 FNEND 

2852 SUB Xyread 
2BsS Shel =e 
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- 260 
Ge oNsies 
«578 
se7s 
288O 
SES 
289 
Bess 
~ 308 
2205 
2312 
2215 
2920 
2325 
e25u 
Fateh oles 
2942 
aos 
2350 
Pagel 
2962 
2555 
PagelTe 
275 
2980 
Zge5 
2390 
Z2a95 
3@22 
3005 
301 
3015 
3020 
3025 
3032 
5035 


5040 


3045 
3052 
3055 
Sb 
3065 
3278 
3075 
3082 
3225 
3032 
S835 
3122 
3185 
ul UG 
225 
3120 
e225 
313@ 
3125 
2140 
3145 
5150 
S155 


INPUT “ENTER FILE NAME" Files 
SEEF 
INFUT “ENTER NUMBER OF X,Y PAIRS" ,N 
ASSIGN @File TO Files 
FOR I=1 TON 

ENTER @File:XxX,Y 

PRINT (X.Y 
NEXT J] 
SUBEND 
SUB Purge 
BEET 
INPUT “ENTER FILE NAME 10 BE DELETED” Files 
PURGE Files 
GOTO 2912 
SUBEND 
SUB Wilson 
COM /Will/ Doa(4),D1a(4),Kma(4),Iact 
COM /Wil2/ Delta,Isat,Nsets,Hod,Cia(3),Alpaa(3) 
DIM EmC4 Smt ss ote 1Ce hea 2e yy al 2g) 
Beer 
INPUT “PLEASE RE-ENTER NAME OF FILE eres 
ASSIGN ®File TO Files 
INPUT “ENTER TUBE €ODE’ ,J]tube 
BEBE 
INE UY = GVEA NAME & OR iy sr ise eK S 
CREATE BDAT Xys,S& 
ASSIGN ®Xy TO Xy$ 
ade oe 
Do=Doa(Itube) 
Di=Dia(Itube) 
Km=Kma( I tube ) 
Ax=PT#01°2/4 |! Cross-sectional area 
Ao=PI*#DoeLl 
Rm=Do#LO0G(Do/D: )/(2#Km) 


Initial values 

Tf=Tsat 

Alpe=.655 

CisSei72 

6=3.651 

Ipeo=vU 

Tend=3 ‘'ICHANGE TO 4, IF FIVE TUBES IN BUNDLE 

IF Tact<1@ THEN 
Ibeg=lact 
Tenc=lact 

ENDS -1F 


FOR IJ=Ibeg TO Iend 
Sx=Q 
Sy=@ 
Sxs=Q 
Sxy=Q 
Jset=2 
ASSIGN @File TO Files 
ENTER @Filé.Dtot, ltube Eorst Od l Does eas sed4 ods 
ENTER @File;:Fm( eo emi Ce inee) 
FOR J=@ TO 8 
T(J )=FNTvsv(Emf( J )) 
NEXT J 
Tvap=(T(Q0)4TC1 )+T(2))/3 
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CY 
tJ 
Wr 
S 


tt OF Ot OI Oy 
tJ $I hU FI ho 


Piro= (2 0+ ia) Z 
IF Isat=@ THEN 
Tsatetliaq 
EloE 
Tsatelvap 
END IF 
Grad=FNGrac( T(I+5S)) 
Delt=ABS( Te’ 1)/(Grad+1Q) ) 
Tavg=I(I+5)+Delt+#.5 


Water/Ethylene Glycol Mixture Properties 
Rhoeg=FNRhoeg! Tavg,Egrat ) 

Nueg=FNNueg!‘ Tavg,Egrat ) 

Mueg=Nueg*+Rhoeg 

Cpeg=FNCpeg/ Tavg,Egrat ) 
Keg=FNKeg(Tavg,Egrat ) 

Preg=Cpeg*Mueg/Keg 


Mdot=FNFmcal(I,T(1+5),Fm(1)) 
Veg=Mdot/(Rhoeg+#Ax ) 

Reeg=Veg*Di /Nueg 
Res=4#Mdot/(PI*Mueg+(D0i-4*Delta)) 
Qdot=Mdot*#Cpeg*Delt 

Qdp=Qdot/fAo 
Lmatd=Delt/LOG((Tsat-T(1+5))/(Tsat-T(1I1+5 )-Delt)) 
Uo=Odp/Lmtd 
Omega=(1+5.484E-3#Preg*.7#(Res/Hod)*1.25)*.5 


R-114 Properties 
Hfg=FNHfg( Tsat) 
Kf=FNK( Tf ) 
Rhof=FNRho( Tf?) 
Muf=FNMu( Tf ) 


Fe(Kf*Z*Rhof*2*GeHfg/(Muf *Do*Qdp ))*.33333 
Ho=Alpa+F 
Two=Tsat-Qdp/Ho 
Tf=Tsat/3+2+Two/3 
Ye(1/Uo-Rm )#F 
X=Do#F/(Keg*+0mega ) 
PRUNI ee OMeors .Omecgd, F=":F. "X=": Xs "Ye": ¥ 
XKa(Jset )=X I SENEEPFIGIENT (“MODIFY LATER ) 
Ya(Jset )=Y 
Sx=Sx+X 
Sy=Sy+yY 
Sxs=Sxs+X#X 
Jset=Jsettl 
IF Jset“Nsets THEN 3]35 
ASSIGN @File TO * 
Slope=(Nsets#Sxy-Sx#Sy )/(Nsets#Sxs5-Sx 2 ) 
Intcpt=(Sy-Slope+*Sx )/Nsets 
Cic=l1/Slope 
Alpac=l/Intcpt 
Cerr=ABS((Ci-Cic)/Cic) 
Aerr=ABS((Alpac-Alpa)/Alpac) 
IF Cerr>. 001 OR Aerr>.@Q@1 THEN 
Alpa=(Alpat+Alpac)+#.§ 
Ci=#(CitCic)+.§ 
PRINT "CIC=";,Cici "ALPA=";Alpa 


83 


fal (2 


6 

46S 
5472 
3475 
2480 
3485 
2490 
3495 
3580 
Sato 
3510 
eis 
5520 
2525 
Sack 
2505 
3540 
3545 
2550 
Bye sie 
SSbie 
S505 
ge 0 
2575 
358 
S555 
3592 
355 
3622 
3605 
3612 
2615 
362@ 
Se25 
3652 
3655 
3642 
3645 
S650 
2655 
3660 
S565 
367@ 
2675 
3682 
3685 
369 
S625 
37Q2 
37@5 
a7 8 
Ses 
et (Ac 
3730 
oG5 
3742 
3745 
3750 
S25 


GOTO 210 
END IF 
Belisle 
Beem 
PRINTEF IS 1 
PR iN Sele Cre. “ALE A= lee 
Crat 1] )=Cic 
Alpaai!l)=Alpac 
PRINTER IS 7@2 
FOR J=@ TO Nsets-1i 

OUTPUT Oxy: Agel Jos yet)? 


NEXT J 
PRINTER IS 1} 
NEXT I 
ASSIGN @xy TO * 
SUBEND 


SUB Nusselt 
COM 7Nus; lime tsat, Gap noe 
Do=.@159 
Ho=12@202@ 
IF JTin=@ THEN 
BEEF 
INPUT “ENTER TSAT AND HEAT FLUX" ,Tsat,Qdp 
END IF 
Hf o=FNHFg( Tsat) 
Two=Tsat-Qdp/Ho 
Tf=Tsat/3+2+Two/3 
K FePNK (TF ) 
Rhof =FNRho( Tf ) 
Muf=FNMUC TF) 
Hoe=.655*( Kf S*Rkhof 2*9.81*hte/i Mut *0e*Gde)) 155554" 
IF ABS((Ho-Hoc )/Hoc )>.@@1 THEN 
Ho=(HotHoc )#.5 
GOTO 3585 
EN Delle 
IF Iin=@ THEN PRINT "“HO="“;Hoc 
SUBEND 
SUB Flot! 
DIM Yaa(4) 
PRINTER IS 7@5 
Idv=} 
BEEF 
INPUT “OK TO USE DEFAULT VALUES (1=DEF=Y,@=2N)",Idv 
IF Idv#l THEN 
Tin=<c 
Xmin=] 
Kmar=5 
Xstep=1 
Ymin=2 
Ymaxe2.@ 
Ystep=.5 
Else 
INPUT “ENTER MINIMUM AND MAXIMUM X-VALUES" ,Xmin,Xma» 
BEEE 


INPUT “ENTER MINIMUM AND MAXIMUM Y-VALUES", Ymin, Yma~ 
BEE 


INPUT “ENTER STEP SIZE FOR X-AXIS",Xstep 
BEEP 


INPUT “ENTER STEP SIZE FOR Y-AXIS", Ystep 
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~J 
~) 


ENS IF 


PRINT “IN: SP1:1IP 2500,1800,8300,680¢@, “ 


PRNGe 50 0, 1OC,OmmbO. TL 2.8." 


Sfx=100/(Xmax-Xmin) 
Sfy=100/( Ymax-Ymin) 
PRINT "PU @,@ PO” 


FOR XKa=Xmin TO Xmax STEP Xstep 


K=(Ka-Xmin)#Sf- 

PRI e res. 5 Os UXT 
NEXT Xa 
PRINT "PA 10@0,8.PU; “ 
PRINT "PU PA 8,8 PD” 


FOR Ya=Ymin TO Ymax STEP Ystep 


Y=(Ya-Ymin)eSfy 
Brunt FOO, aX, YT 
NEXT Ya 
PRINT "PA @,10@ TL @ 2" 


FOR Xa=Xmin TO Xmar» STEP Xstep 


K=(Ka-XmindeSfr 

PRINT PA i:X,°,100; xT" 
NEXT Xa 
PRINT “PA 100,120 FU FA 1208,2 


oe 


FOR Ya=Ymin TO Ymaxz STEP Ystep 


Ye(Ya-Ymin)#Sfy 

Pein orp PA 100," .Y, "XT" 
NEXT Ya 
PRINT "PA 188,100 PU" 
PRINT “PA>@,=2 SR 1.5,2" 


FOR Xa=xXmin TO Xmax STEF Xstep 


X=(Xa-Xmind#Sf» 

PRINT "“PA";X,”,@: ° 

PRINT "CP -2,~-1;LE"ixai"" 
NEXT Xa 
PRINT “PU PA 6,0" 


FOR Ya=Ymin TO Yma» STEP Ystep 


Y=(Ya-Ymin)#Sfy 
PRINT "PA @,";:Y,"" 


PRINT “CP -4,-.25;LB":Yai"" 


NEXT Ya 
BEEP 


INPUT “SELECT MODE (@=HN/H1,1= 


Ism=Ism+1 
IF Idv=l THEN 


IF Ism=1 THEN Ylabel$="HN/H1" 


IF Ism=2 THEN YlabelS$="HN(avg)/H1" 


XlabelS="Tube Number" 
ELSE 
BEEP 


PNEUIe ENTER X-LABEL” ,Xlabel$ 


BEEP 


INPUT “ENTER Y-LABEL" ,Ylabeis 


END SIE 


PRINT “SR 1.5,2;PU PA 50,-1@ CP";-LEN(Xlebel$)/2;:"@iLB"iXlabel$;” 
EA imieere =11,50 CP 0, ~;-LENK Yhabel$)/2*5/6; "01 @,1:LB8"+Ylabeis:"" 


Paints Ce 2.0 DI” 

Okp=1 

BEEP 

INPUT “WANT TO PLOT DATA FROM 


A FILE (1l=DEF=Y,02N)7" ,Okp 
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HN( avg)/H1)",Ism 


4@55 
4062 
4@65 
4272 
4075 
4280 
4085 
4290 
4@95 
4120 
41@5 
4110 
4115 
4120 
4125 
4130 
4135 
4142 
4145 
4152 
4155 
4160 
4165 
4172 
4175S 
4180 
4185 
4132 
4195 
4220 
42@5 
4212 
4215 
4220 

Bes 
4230 
4235 
4240 
4245 
4252 
4255 
4262 
4265 
4272 
4275 
4280 
4285 
4290 
4295 
4320 
43@5 
4310 
4315 
4320 
4325 
4330 
4335 
4340 
4345 
4350 


eisjele 
INPUT "ENTER THE NAME OF THe GA1R Ee reee oii let 
ASSIGN @File TO Ofiles 
BEEP 
INPUT "ENTER THE BEGINNING RUN NUMBER" ,Ma 
BEEP 
INPUT "ENTER THE NUMBER OF X-Y PAIRS STORED" ,Nsets 
Beer 
INPUT “SELECT A SYMBOL FOR THE PLOTTER (1l=+,2=+,3=c.4=0,5=" )" 5S, 
PRINT FU OL 
IF Sy=1 THEN PRINT “SM#" 
IF Sy@#Z THEN PRINT “°SM+" 
IF Sy=3 THEN PRINT “SMc”" 
IF Sy=4 THEN PRINT "SMo" 
IF SyeS THEN PRINT “SM*" 
FOR I=1 TO Nsets 
FOR J=@ TO 3 
ENTER @File:Yaa(J),0 
IF J=@ THEN Ytop=Yaa(Q) 
Yaa( J)=Yaa( J )/Ytop 
NE Xo 
FOR J2@ 10 3 
Ye( J+l1—-Xmin)eSf. 
Y=(Yaa(J)-Ymin)«Sfy 
PRINT “PA",X,Y,"" 
NEXT J 
NEXT I 
Beer 
ASSIGN ®6File TO * 
GOTO 4042 
END IF 
sie lale 
INPUT "LIKE TO PLOT THE NUSSELT RELATION (1=Y,@=N)?7",Oknus 
PRINT "PU; SM" 
IF Oknus=il THEN 
FOR Xa=Xmin TO Xmax STEP ¥stes/S2 
X=(Xa-Xmin)*Sfx 
IF Ism=1 AND Xa>Xmin THEN Ya=Xa’*.75-(¥a-1)°.75 
IF Ism=2 AND Xar>Xmin THEN Ya=Xa7*(-.25) 
IF Xa=Xmin THEN Ya=l 
Y=(Ya-Ymin)#Sfy 
PRIENT = PA Snes ee 
NEXT Xa 
BEEP 
PRINT "PU" 
INPUT “MOVE THE PEN TO LABEL THE NUSSELT LINE VOK 
PRINT “LBNusselt" 
ENO IF 
IF Ism=2 THEN 
Beer 
INPUT "LIKE TO PLOT EXPTL CURVE (1=Y,@=N)" ,Okex 
Nog=@ 
IF Okex=1 THEN 
BEEP 
INPUT “ENTER THE EXPONENT" Ex 
FOR Xea=Xmin TO Xmax STEP Xstep/1@ 
Nq=Nqt+l 
Ya=Xa~(-Ex ) 
X=(Xa-Xmin)#Sfx 
Ye(Ya-Ymin)#Sfy 
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4355 
4568 
4365 
437 
4375 
45 
age 
4390 
asgs 
4400 
44@5 
4412 
441¢ 
4422 
4425 
4432 
4435 
4442 
4445 
4450 
4455 
4462 
4465 
4472 
4475 
4482 
4465 
4490 
4495 
4502 
4505 
4510 
4515 
4520 
bss 
ae 
B35 
4542 
454c 
4552 
45s 
4560 
4565 
4570 
4575 
4580 
45865 
4532 
4535 
4622 
4605 
4612 
4615 
4622 
4625 
4650 
4635 
4640 
4645 


LEUNG mOUSE EO THEN 


PRINTS Pe 


BESeE 


PG uae OM 


PR aNipeioms gt kU 


END IF 


NrExX|] 2a 
BRN PU" 
BEor 


TEU Te MOVE PEN 10 LABEL AND HIT ENTER”, Ob 


PRINT “(Bs=O" 
PRINT “PRe= 10. 
PRINT “LB "clxe 
GOTO 4302 
END IF 

END 1F 

GOTO 4512 

BBEP 


iieote LlRee?TO PLOT KERN RELATIONSHIP (1=#Y,Q=N)?" Yes 


IF Yes=1 THEN 
FOR Xa=Xmin TO Xmax 
Ya=Xa"(-1/6) 
X=(Xa-Xmin)eSf.- 
Y=(Ya-Ymin)¢Sfy 


oN Al ats 


NEXT Xa 
PRINT "PU" 
Beer 
PoPUt “MOVE THE PEN 
PRINT "“LBKern: FU" 
END IF 
ERINT. PU PA 0,@° 
BEEP 


STEP Xstep/ce@ 


ea 


VOREABEE CKERN REUATIONSHIP. 20 


(weve chest PLOT EISSENBERG RELATION (1=¥,Q=N)7",Ohe2 


IF Ohei=1 THEN 
FOR Xa=Xmin TO Xmas 
Ya=.6+.424#Xa(- 
Ket Mek min) Sf >. 
Y=(Ya-Ymin)eSfy 
PRINT “PA’ X.Y, 
NEXT Xa 
PRINT "FU" 
BEEF . 
INF URS BeVJE THE PEN 
PRINT "LBEissenberg: 
END IF 
PRINT "PU SPQ” 
SUBEND 
DEF FNPvst(Tc) 
COieseld/ lft 
DIM K(8) 
IF Ift=@ THEN 
BEEP 


STEF Xstep/1Q 


2p) 


=) Oise 


OSeABEL THE EISSENBERG LINE (Cr 
aes 


fe NieeeVoleCOnRELATION NOT AVAILABLE.FOR R-114" 


STOP 
SNES LE 
IF Lit areeEN 


Oligmea tbo) 2o4556,—-25.00025655,-156,1706546,64.23285504 ,-118.396¢6225 
Wages ete (11722,20.397506/76,1.E9,65 


READ K(+*) 


oy 


4652 
4655 
4662 
4665 
4672 
4675 
4682 
4685 
4692 
4695 
4722 
4725 
4715 
4720 
4725 
4730 
4735 
4742 
4745 
4750 
4755 
4762 
4765 
47702 
4775 
4782 
4785 
4792 
4795 
4822 
4825 
4812 
4815S 
4820 
4825 
4830 
4835 
4842 
4845 
4852 
4855 
4660 
4865 
4872 
4687S 
4882 
48eS 
4892 
4895S 
4920 
4925 
4912 
4915 
4922 
4925 
4332 
4935S 
4942 
4945 


Tat etc (Sate e4e7 25 
Sum=Q 
FOR N=Q@ TO 4 
Sum=SumtK(N)#(1-T)7 (N41) 
NEXT N 
Br=Sum/(Te#(14+K(5 )#( 1-T )4K(6 #0 1-T)°2) C1 -T O/C KC 7 eC 1-T 7% 24+K(8)) 
Pr=EXP(Br ) 
P=?22120020#Pr 
END Sie 
IF Ift*2 THEN 
TfeTc#!],.8+32+459.6 
P=10*(33.20655-43350.98/Tf-3.2635*LGT( Tf )+2.Q0539E-3+T f ) 
END IF 
IF Ift=3 THEN 
A#9.394685-3066.1/(1Tc+273.15) 
P=133.32#10°A 
END iF 
RETURN P 
FNEND 
DEF FNHfg(T) 
COM Ela ari ¢ 
IF Ift=@ THEN 
Tf=T#1.8+32 
Hfg=6.1451958E+1-Tfe(6.9S1Q079E-2+T fell . 3988688E-4+1 .96Q7843E-7eTF) ) 
HfO=Hfg*#2326 
END IF 
IF Ift=1 THEN 
Hf9=2477200-2450#(T-12) 
END TIE 
IF lft=Z THEN 
Tre l# il. St+Sz | 
Hfgo=7.0SS78S7E+1-Tfe(4.8380S2E-2+1.2619048E-4+Tf ) 
Hfg=Hfg#2326 
END Se 
IF Ift=3 THEN 
TRelt+2/ S15 
Hfg=1.35S264E+6-Tke#(6.358Z55E+2+ 1b *. 747462) 
END IF 
RETURN Hfg 
FNEND 
DEF FNMu(T) 
COM /Fld/ Ift 
IF Ift=@ THEN 
ised & Lo 
MumEXP(-4.4636+1011.47/Tk )*1 .E=-3 
END IF 
IF Ift=1 THEN 
A=247.8/(7T+133.15) 
Mu=2.4E-5#107A 
END IF 
IF Ift=Z THEN 
Mu=8.9629819E-4-T#(1.1@94609E-S-T#S.S66829E-8 ) 
END IF 
IF Ift=3 THEN 
Temisv¢{+273<.15) 
MueEXP(-11.0179+TKk8( 1. 744E4+3-Tke(2.80335E+S-Tkel .12661E+8))) 
ENO@ Ie 
RETURN Mu 
FNEND 
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4952 
4555 
4962 
4965 
4970 
4975 
43980 
4985 
43992 
4995 
52002 
5805 
52102 
5015 
5020 
5025 
5030 
5035 
240 
5045 
5050 
SuSo 
5262 
S5@65 
5072 
Su7o 
5082 
5085 
52030 
5095 
5100 
ro) ales 
=) 
SS 
5120 
SkzZS 
Soe 
=) ya}S 
5142 
5145 
5152 
S155 
5162 
Sls 
5172 
Spee 
5182 
5185 
5192 
S135 
5200 
5205 
5210 
Sao 
5220 
5225 
5230 
5250 
5240 


DEF FNVvst( Tt) 
Com /Fid/ lft 
IF Ift=Q THEN 
BEEF 
PRINT “VUST CORRELATION NOT AVAILABLE FOR R-114° 
STOP 
END IF 
IF Ift#l THEN 
PeFNPvst(Tt) 
T=T¢+273.15 
X=1500/T 
Frel/(l1+T#1.&-4) 
F2e( 1-EXP(-X))*2.S5*EXP(X )/K*.5 
B=, 0015#F1-.000942*#F2-. 2004882 +x 
K=2#P/(461.524T) 
Vai 1t(1+2*#Bek )°.S5)/K 
END IF 
IF Ift=2 THEN 
Cre reel. ot52 
Y213.955357-Tfe( .16127262-TF*#S.1726190E-4 ) 
VeV/16.018 
END IF 
IF Ift=3 THEN 
Tk=Gt+2z7s.1S 
P=FNPvst( Tt) 
V2133.95*Tk/P 
END IF 
RETURN V 
FNEND 
DEF FNCp(T) 
COM 7Eig/ 1ft 
IF Ift=@ THEN 
Tezit274e0S 
Co=.40118+Tk#(1.6SQ007E-3+Tk#(1.51494E-6-Tk*6.67653E-10)) 
END IF 
IF Ift=l1 THEN 
Cp#4.211208S8-T#(2.26826E-3-1T9(4.,42361E-S+2.71428E-7*T )) 
ENO IF 
IF Ift=2 THEN 
Co#9.2507273E-1+7T*#(9.3400433E-4+1.7207792E-6eT ) 
ENO IF 
IF Ift]=3 THEN 
Tk®T+273.15 
Co=4.1868+(1.6884E-2+Tk#(3.35083E-3-Tk#(7.224E-6-Tk#7.61748E-9))) 
END IF 
RETURN Cp*1000 
FNEND 
DEF FNRho(T) 
Bom /Eid/ lft 
IF Ift=Q THEN 
Mee 1+275.45 
X£1-(1.8#Tk/753.9S) 
Ro=36.32+61.1464148X*(1/3)+16.41801S5S*xX+17.476838#xX* .St1.119828#x*c 
Ro=Ro/ .@62428 
END IF 
IF Ift=1 THEN 
Ro=999.52946+T+#( .Q@1269-T#(5.482513E-3-T#1.234147E-S)) 
ENO IF 
IF Ift=#2 THEN 
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5245 
5.50 
S255 
S260 
5269 
5272 
5255 
Seeu 
S2eo 
5290 
5205 
5300 
5305 
5310 
e215 
5320 
5525 
5550 
5339 
5342 
5345 
S350 
Sa55 
So50 
Bobo 
eno 
5382 
Saco 
5aa0 
Sogs 
5400 
54@5 
5412 
5415 
5422 
54725 
5430 
S455 
5442 
S445 
5450 
5455 
546¢ 
5465 
5472 
S4a¢ 
5482 
5485 
5492 
5495 
5500 
5505 
5512 
5515 
5520 
5515 
So28 
Sooo 
$542 


Ro=l .GcO7479Ess-1 eC. . 2106 546+i* 2 oe cee 


END =I 


IF Ift=3 THEN 


{ke lee 7S: 


PSS 2oere 


ViwS. 24848E-4474 (6.2 796E-7+1h at S244 de 0+) 4S 7057 e ee 


Ro=1/Vf 
END IF 
RETURN Ro 
FNEND 
DEF FNPr(T) 


PreFNCp(T)*FNMu(T)/FNK(T ) 


RETURN Pr 
FNEND 
COEF ENKGT? 


COM 7Fid7 Tix 
IF Ift#@ THEN K#=.071-.00025i4T 
IF Jft=1 THEN 


MeCT+273. 


[S427 oelS 


w=. 922474+X"(2.8395-X%8( 1]. 8007-4. 52577-07544 ee 


ENGR Le 


TE lit =a 
Kee, 2095238E-2-)#( 2 22142866 449 *2) 30095-2428) 


ENG Se 


Tpit tao (HEN 


Tk=14275.. 


END IF 
RETURN k 
FNEND 

DEE SENET CT } 


LS 


COM ie tealbal ae at 
IF Ift=@ THEN 


BEEF 


PRINT “HF CORRELATION NOT FOR R-114" 


Siar 
END IF 


IF Title hen 
Hr=T*(4.203849-1*(5.881326-4-1 «4. SSiS6Gsi7e-6) ) 


END IF 


feelfti=2 THEN 
rate lae+32 
Hf=6.2078571+Tfe(. 194678S7+Ths] . S2C14286E-4 ) 
pliant 27 Qeyras 


END IF 


lf Jit=S THEN 


Hf=250 | 
ENB @1e 


O° Be VERTETES 


RETURN Hf #1@0@ 


FNENO 
SUB ElotzZ 


COM /Dri/ Star,Sym,Icon 

COM Vela? sift 

OItine Cs svat? 2 loas) 

DATA @.0158,20.0158,0.0158,2.0158 


READ Doa(#) 
Fw=1] 
PRINTER IS 1 
BEEP 
PRINT USING 
PRINT USING 


"4x," "Select Osetion *-Yueimies: ..° 
"6X,°"O Use default values“"” 
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mnm mtn 
ioow ow goteu 
-jry nonin & 
nse nana wm 


SS 


PRUNE MUSING 6x." 
INPUT Ord 

Beer 

INPUT "ENTER TUBE 
Do=Dca(Itube) 
Iht=2 

Seer 


PRINT USING 
PRINT USING 
PRINT USING 
PRINT USING 


maAX 
"BX"? 
mEX 
"BX" 


Sess 


new values. — 


Coda code 


Select, option: 

Qh versus q''” 

1 q versus Delta-T”""” 

< h versus Delta-T (default)""” 


INPUT Iht 
PRINTER IS 7@5 
IF Okd=@ THEN 
IF Iht2=@ THEN 
Ymin2Q 
Yma»x=62 
Ystep=10 
Xmin®.2 
Xmax=1.4 
Xstep=.2 
END IF 
IF Iht2l THEN 
Xmin=2 
Ymin=Q 
Ymax=.S 
Xmax#15 
Xstep=3 
Ystep=.1 
END a 
IF Int=2 
Xmin#®U 
Ymin=Q 
Xmax=52 
Ymax=6 
Xstecp=10 
Ystep=l1 
END IF 
END, or 
IF OkG=1 TREN 
BEE 
INFUT “ENTER 
BEEP 


THEN 


PNEUT eS ENTER 
BEER 
INPUT “ENTER 
Beer 
INFUT “ENTER 
END IF 
BEEF 


IAXIS DEFAULT VALUES 
EChavs. a) 


rCovs t ) 


nevis: 4) 


MINIMUM AND MAXIMUM X-VALUES" ,Xmin,xma~ 
MINIMUM AND MAXIMUM Y-VALUES” , Ymin, Ymax 
emer clzZe FOR X-Axls ,XSteo 


STEP SIZE FOR Y-AXIS" ,Ystep 


PRINT "IN; SP1;1IP 230@,18@@,8302,6800: " 
PRIN? SC (Oe GOR0. 100; Th 220: ” 
Sf»=100/(Xmax-Xmin) 


Sfy=100/( Ymax-Ymin) 


BEEr 
Icg=@ 


iC eleeeTOebYaRASomenGe (]=Y C=N=SODEFAULT )7” Ico 


iF lee eiaeN 6175 


BRING =Pueo,0 PD” 


FOR Xa=Xmin TO Xma- 


STEP Xstep 


ol 


Se4s 
Se8Sa 
Sess 
S860 
S865 
5570 
S87S 
S860 
S885 
S890 
589s 
5900 
S3es 
S910 
S915 
S920 
$25 
5930 
S935 
5940 
S945 
S950 
S9ss 
S960 
SS6S 
5970 
S975 
S980 
S985 
5990 
S995 
6000 
6005 
6010 
6015 
6020 
6025 
6035 
6042 
6045 
6250 
BOSS 
6060 | 
66S 
5270 
6075 
6Q80 
608s 
6030 
6035 
6100 
6105 
kW/m” 
6110 
6115 
6120 
R 1.0,0. 
6125 
6130 


X=(Xa-Amin)*Sf> 
PRINT « PAT tx pee ek ore 
NEXT Xa 
PRINT “FA 100,0.FU, ” 
PRINT “PU PA 0,90 PD" 
FOR Ya=Ymin TO Ymax STEF Ystep 
Yet ya-Ymin?* st y 
RU oleae soy cms evel 
NEXT Ya 
PRINT = PASO, 1@0-1L ao 2° 
FOR Xa=Xmin TO Xmanx STEP Xstep 
X=(XKa-Xmin)*Sfx 
PRINT “PA SX 100-- ht. 
NEXT Xa 
PRINT “PA 100,100 PU PA 100,20 PD" 
FOR Ya=Ymin TO Ymax STEP Ystep 
Y=(Ya-Ymin)#Sfy 
PRINT “PD PA 100,",Y,"YT" 
NEXT Ya 
PRINT —PA 100,100 Fu 
PRINT SPA @.—2 Sk 1. Siz | 
FOR Xa=Xmin TO Xmax STEP Xstep 
X=(Xa-Xmin)#Sfx 
PRINT "PA"; X,",@;" 


IF Xa<i AND Xa<>@ THEN PRINT “CP -1.5,-1;L80:PR =1G eee oe 


IF Xa=@ THEN PRINT “CP -.5,-1:LBQ" 

Xin=@ 

IF Xa MOD 1=@ THEN Xin] 

[fF Xeo=1@ THEN PRINT “CP —=2.—) ee. va ee 


IF Xa>l1 AND Xa<l1@ AND Xin=l THEN PRINT "CP -1.25,-1;LB"iXeai"" 


IF Xa>l AND Xin2@. THEN PRINT “CP -2,-1;LB"ixai"" 
IF Xael THEN PRINT "CP -1,-1:LB1.0" 
NEXT Xe 
PRINT “PU PA @,@" 
FOR Ya=Ymin TO Ymax STEP Ystep 
Y=(Ya-Ymin)eSfy 
PRINT "PA @,";:Y,"" 


IF Ya<i AND Ya<o@ THEN PRINT "CP =4@ — 25° C80°PR —2 0706 3) ene 


IF Ya=@ THEN PRINT “CP -2,-.25.LBQ" 
IF Ya>1 AND InWt<2 THEN PRINT “CP -S,-.25,.LB"iYai"” 
IF Ya>3S AND Int=2 THEN PRINT "CP -4,-.2S5:LE"iYai"" 
IF Ya>@ AND Ye<1@ AND Int]=Z2 THEN PRIN]  €Po-2.-o25 sae 
IF Yael THEN PRINT “CP -4,-.25;1B1.0° 
NEXT Ya 
IF O}d=1 THEN 
BEEP 
[INPUT “ENTER X-LABEL” xiabels 
BEEF 
INPUT “ENTER Y-LABEL",Ylabel$ 
ENO Ie 
IF IhWt<>1 THEN 


PRINT "SR 1.5,2;PU PA =-12,35:D01 @,1;LSn 7 PRele Oo, Ceo PR) Slt boee 


PRINT "PR -1,0.5:SR 1,1.S:082:S5R 1.5,2;PR JSieo ee hee oer ee 


ELSE 


PRINT “PA -12,39;01 @,1;:LBa/(MW/m: SR 1oinS: PRO-10.>c.US2; Shon ee 


Este)” 
ENO IF 
IF Iht=@ THEN 


SV 


eee 


’ co) 
S ? 
o ¢ ;— 


6142 
6145 

AG's 
6150 
2S 
6160 
BLES 
5172 
57S 
6182 
Beles 
61392 
e135 
6200 
6205S 
B21 
6215S 
6220 
BieEZS 
6238 
BcoS 
6242 
6245S 
6250 
Bz 55 
6260 
6265S 
6270 
BZ/9 
6280 
B28S 
6292 
e235 
6302 
6305S 
6312 
1S 
5220 
B25 
E33@ 
ees S 
DIO40 
6345 
6350 
= tsps 
6360 
Be65 
6372 
6575 
6380 
6385S 
6390 
BoegS 
6402 
6485 
6412 
6415S 


Peni Shele Sr wr re 40,-l10,D1l;bBbq/(MW/m PR @.5,1).5R 1,1.5;LB2: SF 
fees, 17.6)" 
ELSE 
Sr ee ech ipa o lanes, Rootes lsrm —2,4,5 PD PR 
PUP R oma: ” 
PRINT "LBwo:PR .S,1;LB)/K" 
END IF 
PRINT “CP 0,0 Ole 
Xlg=l.E+6 
Xug=-1.E+6 
Xal=S2 
Yal=95S 
Nrun=0 
BEEP 
INPUT “WANT TO PLOT DATA FROM A FILE (1#Y,0=N)?7" Ok 
X1l=1.E+6 
Xol==)].E+6 
Okp = 
IF Ok=1 THEN 
BEEP 
INPUT “ENTER THE NAME OF THE PLOT DATA FILE",D_files 
ASSIGN @File TO D_files 
IF Icomb<>@ THEN &26S 
Sx=Q 
Sy=Q 
Sx2=2 
Sxy=0 
Md=! 
BEEF 
INPUT “ENTER THE BEGINNING RUN NUMBER (DEF#=1)",Md 
Noairs®9 
BEEP 
INPUT “ENTER THE NUMBER OF X-Y PAIRS STORED (DEF#=9)°" ,Npairs 
Nrun=NruntNpairs 
PRINTER IS 1 
BEEP 
PRINT USING "4X,°°Select a symbol: "°" 
PRINT USING “4X,°° 1 Star CuaP use some = 
PRINGeUSUNS §4x%,°" 3 Circle 4 Souare’*” 
PRINT USING "4X,"" S Rombus"““ 
PRINT USINGe 4%, 6 Rioght—-side-up triangle ”” 
PRINT USING "4X,"" 7 Up-side-down triangle"”" 
INPUT Sym 
pest il et 
INPUT "ENTER TUBE NUMBER FOR PLOTTING (@=TOP,1=SECOND,...)”,Itube 
PRINTER IS 7@S 
IF Sym=1 THEN PRINT "SMe" 
IF Sym=2 THEN PRINT "SM+" 
IF Sym=3 THEN PRINT "SMo" 
IF Md>1 THEN 
FOR I=} TO (Md-1) 
ENTER @Fileixya(+) 
NEXT I 
END IF 
FOR I=] TO Npairs 
ENTER @FileiXya(*) 
Ya=Xya(Itube#2 ) 
Xa=Xyae(ltube#2+1 ) 
Yc=LOG( Xa) 


Se: 


mManomam 
bP ph ah se 
& B OF OF hot) 
in «& Ul gd Ul FD 


6452 
6455 
6460 
6465 
6470 
6475 
E480 
6485 
6492 
6495 
ESee 
6505 
E510 
BelS 
6520 
sya 
6550 
655s 
£540 
6545 
Go5U 
GS55 
6562 
5565 
65702 
Biss 
6582 
6535 
6530 
6Se5 
560e@ 
66@5 
6616 
66iS 
BEL2 
6625 
6620 
6635 
6642 
6645 
6652 
6655 
5660 
6665 
bE75 
5682 
6685 
6690 
Goae 
6782 
6705 
67g 
bv is 


AC=LOO(Xa/ Ya) 
CxreSxtrXe 
Sy=Sy ee 
Sx cC@ov ee ec L 
Sry=Ssytacryc 
IF Inht=@ THEN 
Xt=Ya 
Ya=Ya/xXe 
Ka=At 
TF Xa/].©+6>xXu1] THEN XulexAa/i-e+6 
IF Xa/1.E+6<x%11 THEN X11]=xea/1.E+6 
Eno 
IF Jht=@ THEN 
X=(Xa* | 26-6 -i mine St 
Y=(Yael .E-3-Ymin)#Sfy 
END IF 
TF Iht=1 THEN 
X=(¥a-Xmin)*Sf- 
We Yee eo Mi he Sy 
END IF 
Le ae ae 
Metta, YaaXmi ty )* ot = 
Y=(C¥a* lob -5- ymin 2esty 
ENDS LE 
IF Y>10@ OR Y<@ THEN 6525 
IF Sym>3 THEN PRINT Sm 
IF Sym: 4 THEND PRINT “Sheol 34 
PRINT =P Aer ka 
IF SymoS THENSERINY “Skies ioe 


IF Syme THEN PRINT “UG2,4,S900,5e.-4,0,0,0.4,05) 


IF Sym=S THEN PRINT “UC3,0,99,-3,-6,-3,6,3,.6,3.-63° 


IF Sym=6 THEN PRINT 


"UGO, 5.0. oo ome es OU ome 


[esymey THEN PRINT )°UCO(-5.5,99,55,0 040. = om cum 


NEXT 1] 
BEEP 


INEUT “WANT TO LABEL. Clay, G=N 7 io] 


LF ffel=1 THEN 
IF Sym>3 THEN PRINT 
IF Sym<4 THEN PRINT 
PRINT "PA", Xal,Yal, 
IF Syn So He Nee aN 
IF Sym=4 THEN PRINT 
IF Sym=5 THEN PRINT 
IF Sym=6 THEN PRINT 
IF Sym=7 THEN PRINT 
PRINT "SM" 
IF Sym<4 THEN PRINT 
PRN) “PR Zyl. or 
Yal=Yal-5 
BEEF 
IF Tas=1 THEN 
BEEP 
INPUT “ENTER THE 
PRINT "PR 2,@:SR 
GOTO 6665 
END IF 
END IF 
BEEP 
INPUT "WANT TO COMBINE 


"em" 
USR Ol i84 ciao 


ajow 


“SO Ra doe selena 

"UGiy ee oonO eto — 4.4050 Ona. Oe 
"UCS,; 05.9055, -O1-o Ojo weno. aoCms 
"UCO , 5.3799) eae eon eos Com 
"UCO,=5.5,989 S76" Oy Oye Suen. 


"PR 2,0" 
150,1.8;LB8° : Oa lesa 


STRING" ,Label$ 
1,0,1.8° CBs babels 3 


ANOTHER FILE? (1=Y,@=N)" ,Icomb 


6720 
6725 
6730 
6725 
6740 
6745 
6752 
Biss 
6760 
6765 
6772 
5775 
6780 
6785 
679 
6735 
6822 
6805 
6810 
6815 
6820 
6825 
6830 
6835 
6840 
6845 
6852 
6855 
6862 
6865 
6872 
6875 
6882 
6885 
6892 
6895 
6922 
69@5 
63912 
6955 
4 
6925 
6932 
G91S5 
6940 
6945 
6952 
63555 
6962 
S965 
6972 
6975 
6982 
6985 
6992 
6335 
(000 
7085 
7@10 


ASSIGN ®@File TO * 
X11=5 
Xul=45 
IF Icomb<>@ THEN 6222 
BEEP 
INPUT “WANT TO PLOT A LEAST-SQUARES LINE (1=Y,Q@=N)",I]s 
IF Ilse! THEN 
BEEP 
INPUT "SELECT EXPONENT: Q=COMPUTE, 120.75)" ,lIexp 
BEEP 
INPUT “SELECT CURVE TYPE (@=SOLID,1=DASHED)" ,I1t 
Tit=lI1t+! 
PRINT "SM" 
IF Iexp=Q@ THEN 
Bb=(Nrun#*Sxy-Sy#5x )/(Nrunt*5x2-5x%2 ) 
ELSE 
Bb=.75 
END IF 
Ae=(Sy-Bb*Sx )/Nrun 
Aa=EXP( Aa) 
PANTER 15 ft 
PRINT USING “10x, “awe "",Z.4DE";Aa 
PRINT USING "10@X,""n ' 2. 40DE ; Bb 
PRINTER IS 705 
In=@ 
IF Iht=@ THEN Xxstep=Xstep/4Q@ 
IF Iht>@ THEN Xxstep=Xstep/1Q 
FOR Xa=Xll1 TO Xul STEP Xxstep 

IF Xa>.99*Xmax THEN 6995 

IF Iht=1 THEN Ya=Aa*xXa*Bb 

IF Iht=@ THEN Ya=Aa*(1/Bb)*#(Xael.E+6 )*(( Bb-1)/Bb) 

IF Int=2 THEN Ya=Aa*Xa*(Bb-1 ) 

IF Iht#=@ THEN 
Y=(Yael.,E-3-Ymin)*#Sfy 
X#(Xa-Xmin)#Sfx 

END IF 

IF Iht=1 THEN 
Y=(Yael.E-6-Ymin)*Sfy 
X=(Xa-Xmin)*Sfx 

END IF 

IF Iht=2 THEN 
Ye(Yael .E-3-Ymin)¢Sfy 
X=(Xa-Xmin)*Sfx 

END IF 

IF Y<@ THEN Y#@ 

IF Y>1@0 THEN 6992 

IF Ilt=1 THEN 
PRIN? “PA X,Y, PO" 

ELSE 
IneIn+l1 
Ir=In MOD Ilt 
if Ir=1 TRENDPRINI "PA" .X,Y,°PD° 
IF Ir*@ THEN PRINT “PA",X,Y, "PU" 

END IF 

NEXT Xa 
PRINT "PU" 
ENU {tr 
Icomb#@ 
Sng) sists 


V5 


7015 
7020 
(sires 
7030 
7035 
7042 
7045 
7050 
7055 
7060 
7@65 
7070 
no75 
7080 
7085 
7090 
¢o3s 
71QQ 
7105 
ees’ 
qi15 
7120 
Wi2s 


7iLse 


i235 
7140 
7145 
R150 
w155 
7160 
Ui 
7170 
ees 
Fee} 
Fass 
7190 
e295 
7200 
7205S 
7210 
el5S 
VE20 
@z225 
(250 
@2o5 
7242 
7245 
7250 
Weo5 
7260 
i255 
7270 
a215 
7280 
232 
wesc 
feds 
7300 
7o05 


END IF 
PRINT °PU SM" 
BEEP 
INPUT “WANT TO PLOT NUSSELT LINE (1=Y,@=N)7° , Inp 
IF Inp=@ THEN 7125 
Beer 
INPUT "ENTER TSAT (DEFAULT=18 DEG C)",Tsat 
Hf g=FNHfg( Tsat) 
X1lleS 
Xul=45 
FOR Xa=Xl1] TO Xul STEP Xstep/5@ 
Tfilm=Tsat-Xa*.5 
KfeFNK(TFilm) 
Rhof=FNRho(Tfilm) 
Muf=FNMu( Tfilm) 
Yat. 7Z284(Kf°S*Rhof*2*9.8l*Hfe/(Muf*Do#xa) )° .25 
X=(Xa-Xmin)*Sfx 
Yet Yar! -€-s-ymind*Ssfy 
PRINT "PA" ,X,Y, "PD" 
NEXT Xa 
PRINT “PU PA @,Q" 
PRINT "PU PA @,@ SPQ" 
SUBEND 
SUB PlotZ 
COM /DOrl/ Star ,Sym. leon 
CORE Vay salt 
Dit cs!) 
Fw=1 
PRINTER IS 1 
BEEP 
PRINT USING “4X,° "Select Cptitonw-) sermtts. 
PRINT USING "G&X,""@ Use default values""" 
PRINT USING "6X,""1 Use new values""" 
INPUT Okd 
PRINTER IS 785 
IF Okd=@ THEN 
Xmin=Q 
Ymin=@ 
Xmax=15 
Ymax=15 
Xstep=3 
Ystep=3 
ELSE 
BEEP 
INPUT "ENTER MINIMUM AND MAXIMUM X-VALUES”,Xmin, Xmax 
BEEF 
INPUT “ENTER MINIMUM AND MAXIMUM Y-VALUES" ,Ymin, Ymax 
BEEP 
INPUT “ENTER STEP SIZE FOR X-AX1IS" ,Xstep 
BEEP 
INPUT "ENTER STEP SIZE FOR Y-AXIS" ,Ystep 
END IF 
Been 
PRINT “IN; SP1;IP 2300,1802,8300,6800; " 
PRINT ~S600,100.0-100 (iie2-O,- 
Sfx#10Q0/( Xmax-Xmin) 
Sf y=100/( Ymax-Ymin) 
BEE 
Icg¥Q 


6 


7310 INPUT “EIKE FO SY=PASSeCAGE (1=Y,CeN=BEFAULT)’” leg 
cS IF Icg=1 THEN 7620 
e2vU PRIN «FU 0,0 PD” 


i FOR Xa=xXmin TO Xmaz STEP Xstep 
7355 PEIN] "PA" sxe 0; XTi” 
q2A0 NEXT Xe 


7345 PRINT FA LOO. Orel.” 
7250 PRINT) sPURPAOONErD © 


7355 FOR Ya=Ymin TO Ymax STEP Ystep 
7360 Y=(Ya-Ymin)*Sfy 
7205 PRUNG eG Oe arty. YX le 


Vai? NEXT Ya 
7375 PRINT "PA @,10@ TL @ 2" 


7380 FOR Xa=Xmin TO Xmax STEP Xstep 
7385S X=(Xar~-Xmin)*#Sfx 
7392 PRINT “PA's Xe" 100, x1” 


7255 NEXT Xa 
7422 FRINT "PA 10@,10@ PU PA 120.0 PD" 


7405S FOR Ya=Ymin TO Ymax STEP Ystep 
7412 Y=Cya-imin Sty 
7415 PRINT "PD PA 100,",¥,"YT" 


7620 NEXT Ya 

ga. S PRINT "FA 100,10@ FU" 

743 ERIN “FAO. c SR 1.5.20" 

7435 FOR Xa=Xmin TO Xmar» STEF xstep 


7440 X=(Xa-Xmind#Sfx 

7445 PRINT “PA";X,",@:° 

7452 IF Xa<l AND Xav>@ THEN PRINT "CP -1.5,-1;LB@;PR -1,0:LB";:ixXai"" 
74S IF Xa=@ THEN PRINT "CP -.5,-1:LB0" 

7462 Xin=Q 

7465 IF Xa MOD 12#@0 THEN Xine=l 

7470 IF Xa>=1Q0 THEN PRINT "CP -2,-1,LB";Xai"” 

7275 IF Xav1 AND Xas1I@ AND Xin#=1 THEN PRINT "CP -1.25,-1,.LB°:xXa,"" 
7482 IF Xa>l AND Xin=@ THEN PRINT "CF -2,-1:LB":iXai"" 

7485 IF Xa=l THEN PPINT "CP -1,-1:LB1.Q" 


7492 = NEXT Xa 
7495 PRINT "PU PA 0,0" 


7520 Tht=2 NOG Ley 

Voe@5 FOR Ya=Ymin TO Ymax STEP Ystep 

1216 Y=(Ya-Ymin)e¢Sfy 

Faas PRINTeORAr@. "3 Y., 

q say IF Iht=@ AND Yav@ THEN PRINT "PR 2,0" 

7525 IF Ya<1 AND Ya<>@ THEN PRINT "CP -4,-.25;LBQ0;PR -2,0:LB";Ya:i"” 
7530 IF Ya=@ THEN PRINT "CF -2,-.25:L5Q" 

oS IF Wa ean int<2 THEN PRINT "CP —-5,-.285;:LB";Ya:"* 

7542 LE asa Oomiat=2 THEN PRINT ~CP —4,-.25;LB": Ya: "" 

7545 [ever Cmeaneera oO aND ITht=2 THEN PRINT “CP -3,-.25;LB":Ya:"" 
7550 | IF Ya=1 THEN PRINT "CP -4,-.25;LB81.0" 


(S55 NEXT Ya 
7560 IF Okd=@ THEN 


7365 XlabelS="X" 

tan Ylabel$S="y" 

Tors BibSr 

758 BEEP 

71385 INPUT “ENTER X-LABEL” ,Xlabel$ 
7590 BEEP 

7335 INPUT “ENTER Y-LABEL”,YlabelS 


7622 END IF 
7605 Poi ohel. Seer UePoeouy 10 0c 1-LeN(Xlabel$)/2:"O:LB°:XlabelS;:”"” 


7618 
7615 
7622 
(B25 
7632 
(652 
7648 
7645 
7550 
(555 


foow 
(66¢ 
7678 
iets 
7680 
1685 
7592 
toes 
776@ 
Tis 
T71E 
7715 
7728 
1725 
7728 
VsS 
774Q 
C145 
775@ 
Tis2 
7762 
(isha 
7772 
tas 
7786 
765 

7739 
7735 
7502 
7805 
7E1O 
Tos 

7820 
Tees 

7650 
teo= 

7842 
7845 

785@ 
tess 

785@ 
7865 

7878 
(37s 

7882 
7885 

7898 
1e32 

790@ 


PRINT “PA -11,5@ CP Q,".-LENCY¥lebelSo7aea76, Ul O) boy epee 


PRINT “CP 0,6 D1” 
Nrun=@2 


BEER 


INPUT “WANT TO PLOT DATA FROM A FILE (1=#Y,0=N)7",Ok 


Okp=@ 

IF Ok=1 THEN 
BEEE 
INPUT 
ASSIGN @File TO D_files 
Ir LCOMD<S>U THEN /b3U 
Sx=@ 

Sy =@ 
Sx2=0 
Sxy=@ 
Md=1] 
BEEe 


“ENTER THE NAME OF THE PLOT DATA FILE” ,D_files 


INPUT "ENTER THE BEGINNING RUN NUMBER (DEF=1 >)" ,Md 


Npairs=9 
BEEe 
INPUT 
Beer 
INPUT 
Nrun=NruntNpairs 


PRINTER -PSet 


“ENTER THE NUNGER OE x 


FAIRS STORED (8EF—2>  INeaiee 


“SELECT TUBE NUMBER (@=TOP,1=SECOND (ee seu. 


SEER 

PRINT USING “4X > “Select te esympo ts. 

PRINT MUSING “4X 4092) SS tec 2 SF luisesagn im 
PRINT TUSING "4X5" ss Cire2 eu Squere’"" 
PRINT USING "4X,"" S Rombus""" 

PRINT USING "4X,"" 6 Right-side-up triangle"”” 
PRINT USING "4X,°"° 7 Up=sicescdoun treancie 2: 
INPUT Sym 

PRIN eR le )7es5 

IF Sym=1 THEN PRINT "SMe" 

IF Sym=2 THEN PRINT “SM+" 

[FP eSym=s THEN PRINT) = Sho. 


Ma=Itube*9 
IF Md-l THEN 
FOR 120 TO (Md-1) 
ENTER @Fileix,Y 
NEYT J 
END IF 
FOR I=1 TO Noairs 
ENTER @File:irxa,Ya 
Sx=Sx+Xe 
Sy=Sy+Ya 
Se Soa 2 ae 
Sxy=Sxyt+XarVYa 
X=(Xa-Xmin)eSfx 
Y=(Ya-Ymin)eSfy 


IF Y>1@@ OR Y<@ THEN 7910 


IF Sym>3 THEN PRINT 
IF Sym<4 THEN PRINT 
PRINT “PA",X,Y,"" 

IF Sym>3 THEN PRINT 
IF Sym=4 THEN PRINT 
IF Sym=S THEN PRINT 
IF Sym=6 THEN PRINT 


mu Sm” 
MOR 1 Oe a4. 


"GR isd ee 

"UC2 4,990.58 .=4 Cyne ace 
"UC3..\0 SS = oreo ye Sib sba oS oe 
"UC On S. S555 54564 moO oa 


9% 


79@5 
7910 
al 
7928 
We © 
7932 
7255 
7940 
7945 
7950 
T2599 
7962 
7365 
7970 
7 elfes) 
7982 
Ga85 
feziche| 
8Q22 
8ee@s 
8212 
8@15 
BAce 
S025 
8230 
8035 
8240 
8245 
8252 
8@S5 
8262 
8@65 
8272 
Be75 
8280 
8@85 
8290 
8295 
812¢ 
€12@5 
8110 
es 
8122 
625 
2132 
@13s 
8142 
6145 
8150 
shUSsS 
8162 
Sles 
8178 
el 76S) 
8182 
Big5 
8192 
8S 
8220 
8205 


Precoyh= 7 (Henman CCCs 5.5. Jose oO. - 5,-6.5.” 
NEXT I] 
Beer 
TNPUT “WANT TO LABELOIl=Y .@=No7 > 2115) 
IF I]lbl=1 THEN 
lee ayn TOeNeE REN a: om 
IF Sym<4 THEN PRINT “SR 1.4,2.4° 
Pein SFA eke beyal...”” 
lewoyM7 Ss TRENT PRINT *SR LiZ2e ies 
IF Sym=4 THEN PRINT “UC2,4,99,2@,-8,-4,0,0,8,4,0:" 
tpeove= shew rhunie UCS, 0, d0.—-o8-b.-o.6,6.60,5,—-6:; " 
[Ewove=Oe linens RINT “UCO,5.3,99.3,-8,-6.2,3.8: " 
Te Sym=7 THEN PRINDUEO, -5.3,99,-3,.8,6,9,-2,-8:’ 
PRINT “SM” 
IF Sym<4 THEN PRINT “PR 2,2" 
PrOnNGeerRec 1 .0;5n0! 0.1.8; LB° 7 O_ files: ~* 
Yal=Yal-5 
INPUT "WANT 70 ADD ANOTHER STRING (1=Y,@=N)?", las 
IF Tas=l1 THEN 
Beer 
TNEUT “ENTGr THE STRING’ ,Labels 
PROUNT  PRaeGnonk leo. lS ,e «labels: >= 
GOTOs7 390 
END IF 
END IF 
Beer 
INPUT "WANT TO COMEINE ANOTHEF FILE? (1=Y,@=N)"“,Icomb 
ASSIGN @File 70 * 
IF Icomb<>@ THEN 7645 
Ils=1 
Sele 
INPUT "WANT TO PLOT A LEAST-SQUARES LINE (]=DEF=YES ,@=NO)", Ilse 
IF I]ls=1 THEN 
BEER 
Phew =seeec) CURVE TYPE <@=S0LID,1=GASHEL?” ,1T1t 
Mae 11 t+4 
PRN eo 
IF Texp=@ THEN 
Boe Neum*Sxy-Soy*Se 07 (Neun*Sac-Ss 2? 
PESe 
Bbh=.75 
END IF 
Aa=(Sy-6b*S» )/Nrun 
PRINTER IS 1 
POR eUSuNGe 10x, “vas= "" M2 eeDE”: Aa 
PRINT UsitNG “10x%.° on © Fi iae SOS) 0: 
PRINTER IS 785 
In=@ 
FOR Xa=Xmin TO Xmax STEP (Xmax-Xmin) 
Ya=Aa+xXarBb 
Y=(Ya-Ymin)#Sfy 
X=(Xa-Xmin)*Sf-x 
IF Y<@ THEN Y=@2 
IF Y>10@ THEN GOTO 8220 
ie t= 1 THEN 
ms UN ho = letra Qe fen Wl Big 
EPse 
In=In+l 
ie=ln MO Tit 
fe e= le iheN PRINT PA’. X,Y,” PD” 


29 


IF fe=0 THEN PRINT ee eee? Lie 
END IF 
NEXT 
PRU Niet Use 
END IF 
Icomb=2 
GOTO 7620 
END IF 
PRINT "PU PA @,@" 
PRINT "PU PA @,@ SPQ" 
SUBEND 


1900 
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